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Abstract 

Theoretical aspects of x-ray standing wave method for investigation of the real structure of crystals 
are considered in this review paper. Starting from the general approach of the secondary radiation yield 
from deformed crystals this theory is applied to different concreate cases. Various models of deformed 
crystals like: bicrystal model, multilayer model, crystals with extended deformation field are considered 
in detailes. Peculiarities of x-ray standing wave behavior in different scattering geometries (Bragg, Laue) 
are analysed in detailes. New possibilities to solve the phase problem with x-ray standing wave method are 
discussed in the review. General theoretical approaches are illustrated with a big number of experimental 
results. 
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1 Introduction 



A new field in the physics of x-ray diffraction has appeared and successfully developed during last 30 years. 
It is based on studying and using x-ray standing waves (XSW) that are formed in a perfect crystal under 
conditions of dynamical diffraction. Apart from general physical interest involving the enormously sharp 
change in the interaction of x-rays with atoms in the crystal and on its surface, this field, as has now become 
clear, is highly promising for analyzing the structure of crystals and its adsorbates at the atomic level. 

Actually a standing wave that has the same period as the crystal lattice is extremely sensitive to the 
slightest deviation of the atomic planes (or individual atoms) from their correct position in the perfect crystal 
(or on its surface). Thus XSW method is particularly useful in its application for structural analysis. For 
this technique, an x-ray interference field (XIF) is produced by the superposition of, typically, two plane 
x-ray waves. In this case we have the following expression for the amplitude of the electric field in the crystal: 



E(r) 



E e ,;k ° r 



ik h r 



(1) 



where ko is an incident wave vector, = ko + h, and h is the reciprocal lattice vector multiplied by 2tt. 
The field intensity is determined by the square of the modulus of the amplitude E(r) and is equal to 



J(r) 



lEnl 2 



1 



IE, 



2^cos(hr + a(0)) 



(2) 



where a(9) is the phase of the ratio E/j/Eq. The spatial position of the planar wave field is determined 
by the phase a{9) between the two (electric) field amplitudes. Generated via Bragg reflection, employing a 
diffraction vector h, the x-ray standing wave exists within the overlap region of the incident and reflected 
x-ray wave (Fig.l) and the phase a{9) and thus the position of the wave field is a function of the angle 
9 measured from exact Bragg angle, varying by half a diffraction plane spacing within the total reflection 
range. Thus, atomic positions can be scanned by the XIF and exactly determined if the yield of the element 
specific photoelectrons or x-ray fluorescence photons is recorded as a function of the glancing angle. 

Structural analysis by the XSW technique represents actually a Fourier analysis but, in contrast to 
diffraction techniques, the atomic distribution of an elemental sublattice is sampled. The two important 
parameters which are determined by an XSW measurement are called coherent fraction (F^) and coherent 
position (P]}) and represent the h-th amplitude and phase, respectively, of the Fourier decomposition of 
the distribution of atoms under consideration. The XSW method is particularly powerful for the analysis 
of the structure of adsorbates on crystalline substrates since the position of the adsorbate atom within 
the surface unit cell can be determined with high accuracy for low adsorbate coverages. In case several 
elements are present on the surface, F c ' 1 and P£ can be obtained for each elemental sublattice within one 
XSW measurement. 

An effect involving the existence of a standing wave and the variation of the total field at the atoms of 
the crystal lattice has been known for a long time (see for example (l], |^, || ) . However in the conditions of a 
classical x-ray diffraction experiment, when the intensity of the reflected and transmitted waves is measured 
separately it is manifested very weakly. This is mainly due to the fact, that the cross sections of the inelastic 
scattering channels are considerably smaller, than the cross section of elastic scattering The standing 
wave in the crystal reveals itself in a traditional x-ray diffraction experiment only in the form of an anomalous 
angular dependence of the absorption (the anomalous transmission effect in the Laue case as discovered by 
Borrmann §|) and also the weak asymmetry of the reflectivity curve in the Bragg case jjj 0, E§. 

Batterman || was the first who made an attempt to see the standing wave and its behavior by measuring 
the GeK a fluorescence emitted by crystal atoms. Despite expectations, the measured curve has the angular 
dependence similar to the inverted reflectivity curve of the x-rays. The structure of the wave field manifested 
itself very weakly only at the edges of the total reflection region. It was soon understood , that this behavior 
was due to the fact, that the depth of yield L y i of the fluorescence radiation exceeds by far the penetration 
depth of x-rays in the crystal. This penetration depth while the dynamical diffraction of x-rays is of the 
order of extinction length L ex . As a result all the radiation absorbed in the crystal gives rise to a fluorescence 
signal. Its amount, following the law of conservation of energy is equal to 1 — Pr(9), where Pr{9) is the 
reflectivity curve. The secondary radiation (SR) yield is proportional to the wave field intensity (^) at the 
atoms only if the condition (L y i « L ex ) is fulfilled. Later there were proposed methods for revealing the 
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structure of the wave field by measuring the fluorescence yield from the impurity atoms introduced in the 
lattice of the crystal matrix at a very small depth Q or measuring the fluorescence signal at grazing exit 
angles 0. Evidently the condition L y i < L ex is satisfied in this cases. Moreover, for a monolayer of atoms 
absorbed on the surface of crystal this condition for the escape depth of the SR is surely satisfied. 

The above mentioned problem does not exist in measurements of the photoclcctron emission since elec- 
trons escape from a thin subsurface layer with a thickness of fractions of a micrometer. Already in the first 
works on the measurement of the photoclcctron emission, carried out in the former Soviet Union in the 
early 70 's g |k}, the dispersion like angular dependence corresponding to the behavior of an x-ray standing 
wave (^) was observed. What was understood from the very beginning that this angular dependence of the 
photoelectron yield curve contain essential information about the structure of the surface layers. Later this 
field of research was developing intensively in several scientific centers of the former Soviet Union (see for 
details review paper Q and a book |^2| ). 

Already in the middle sixties first attempts to measure different secondary processes were made. For 
example measurements of thermal diffuse and Compton scattering while the existence of the standing wave 
in the crystal were reported (k| [ll| (see also later experimental paper jl6| ). The angular dependence of 
the photoelectric current in the silicon crystal with p — n junction while the dynamical scattering of x-rays 
was measured Jf7], [T^] . 

During last decade, due to the availability of the synchrotron radiation facilities of last generation (ESRF, 
APS, Spring-8) XSW method has become a useful and even in some cases a routine tool for investigating the 
surface of the crystals and the structure of the adsorbates. Most of the results obtained up to the beginning 
of 90-th, especially applications of XSW technique to surface analysis were summarized in a review paper 

, an overview of the method was also given in a number of papers |2(| [|| . However one of the 

important field of applications of the XSW method for investigation of the structure of the real crystals 
(containing different type of defects, implanted crystals, epilayers on the surface of the perfect crystals, 
geterostructures etc.) has not been reviewed up to now. Previous review on this subject was written 
nearly fifteen years ago (see also the book and a big number of new results are not summarized until 
now. At the same time still it is a big interest to the foundations of the theory of XSW in real crystals (see for 
example recent paper |24[]). In our work we are planning to fill this gap. Theoretical approach is illustrated 
by the experimental results obtained in the Laboratory of Coherent Optics and Synchrotron Radiation of 
the Institute of Crystallography RAS. Due to a limited size of this manuscript we have no possibility to give 
a detailed consideration of all results obtained in the field of XSW method in different research centers all 
over the world. Some of them are just mentioned or even not mentioned, but this, surely, does not mean 
that they are not relevant to the subject. This can be a special subject of another review paper. 

If XSW method in perfect crystals is based on the dynamical theory of x-ray diffraction (see for e.g. 
books and reviews |], ||, ||, [2^, 26 1) for the description of the fields and the yield of the secondary radiation 



from the real crystals it is most effective to use Takagi-Taupin theory [[27], |28|, |29[ of the propagation of 
x-rays in the deformed crystals. For convenience of the reader we start Chapter II with formulation of the 
main results of this theory that will be used in the following parts of the work. In the end of the same 
Chapter for the same reason we give the main relationships for the description of the real and imaginary 
part of the susceptibilities in crystals in x-ray wavelength region. For a recent review of the dynamical 
theory of x-ray diffraction in a perfect and deformed crystals see also Q, pi] . Chapter III gives a general 
mathematical formalism for calculating the secondary radiation yield in a real crystal. This chapter is based 
on the results of the paper J32] | and represents the theoretical foundation for the remainder of this review 
paper. Next Chapter IV is devoted to the theory of XSW in the case of perfect crystals. Peculiarities of the 
secondary radiation yield with the big and small depth of yield of the secondary radiation are discussed, in 
the next subsection fluorescence and photoemission yield from the crystals containing different type of atoms 
is analysed and in the end of the chapter crystals with amorphous surface layer are discussed. Chapter V 
is devoted to very important and often realized case of deformed crystal, that can be approximated in the 
frame of bicrystal model. In the first subsection theory of the secondary radiation yield from such a model 
crystal in the Bragg geometry is presented and in the following subsection it is illustrated by a numerous 
examples. In the next Chapter VI Laue geometry is considered, peculiarities of x-ray standing wave behavior 
in this geometry are discussed and illustrated by examples. In Chapter VII the bicrystal model is generalized 
to the case of a multilayer model of the deformed layer and a secondary radiation yield from such system 
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is analyzed theoretically and with its applications to the study of implanted crystals. Secondary radiation 
yield from the crystals with extended deformation yield are considered in Chapter VIII. There is given a 
detailed description of the wave fields in the case of the crystals with the uniform strain gradient, which 
includes the case of bent crystals and as a special case vibrating crystals. Next Chapter IX describes one 
of the important applications of the XSW analysis: the possibility to solve a phase problem while x-ray 
scattering from deformed crystal. This approach opens the possibility to determine uniquely the structure 
of the surface layer directly from the scattering experiment. Last Section X presents a summary and an 
outlook for the future applications of the XSW method in real crystals. 
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2 X-ray dynamical diffraction in real crystals 



2.1 Takagi-Taupin equations 

Directly from the Maxwell's equations for the electric field vector E(r, oj) (oj is the frequency of the incident 
wave) inside a crystal we can obtain the following wave propagation equation, 

(A + k 2 )E(r,w) - graddivE(r,oj) = -k 2 — j(r,w), (3) 

OJ 

where k = |k| = oj/c is the magnitude of the wave vector (c is the velocity of light), j(r, oj) is the current 
density induced by the electromagnetic field. This current in the case of linear electromagnetic wave theory 
is, in fact, a linear function of E(r, oj), 

ji(r,ui) = J dr , a ik (r,r',oj)E k (r',oj), (4) 

where ai k (r, r', uj) is the nonlocal tensor of the conductivity of the crystal. In general case equation (^) 
describes non-local coupling between j(r, oj) and E(r, oj). It takes into account all possible interactions (such 
as elastic Thompson scattering, photoelectric absorption, Compton scattering and an inelastic scattering 
on thermal phonons) between the electromagnetic wave and the crystal [ j33[ . The main contribution to 
<7ifc(r, r', w) is connected with elastic Thompson scattering and has a strictly local character (the same is 
valid for the main inelastic process, that is photoelectron absorption in the dipolc approximation) so we can 
present the tensor of the conductivity in the following way: 

o- 4fc (r, r', oj) = cr(r, oj)5 lk 5{r - r'), (5) 

where 8i k is the Kroneker symbol and <5(r — r') is the Dirac 6— function. 

For further consideration, if it is not specially noted, we will assume local coupling (|J). According to (||), 
the right hand side of equation (||) takes the form, 

j(r,«) = x(r,w)E(r,«), (6) 

OJ 

where x(r,cj) = (4m/u>)a(r, u>) is the crystal susceptibility, related with the permittivity e(r,oj) of the crystal 
by usual equation: e(r,oj) =1 + x(r,w) []. 

In a perfect (ideal) crystal susceptibility x(r) is a periodic function with the period of the crystal lattice 
x(r) =x( r + a )i a is the translation vector. It can be therefore expanded as a Fourier series, 

X (ld) (r)=E^ d)cx P(* hr )' ( ? ) 

h 

where h =27rH, H is the reciprocal lattice vector. 

We shall assume now, that some part of a crystal lattice (in most of applications it is a thin surface part 
of the crystal (see Fig. 2)) is weakly deformed due to epitaxial growth, implantation or some other type of 
deformation or defects. It is convenient to describe this weak deformation field of a crystal lattice by two 
functions. The first one is the deformation vector u(r), which determines the displacements of atoms in a 
crystal from the position of perfect lattice and the second one is the static Debye- Waller factor e~ w ^ which 
takes into account the random displacements of the atoms from the equilibrium positions in the h direction. 

In the case of weak deformations, that means that relative displacements are small on interatomic dis- 
tances, 



«1, (8) 

dx k 

the susceptibility of the crystal x( r ) is defined from that of a perfect one according to relation [ p8[ , 

X(r)=x (id) (r-u(r)). (9) 



1 Further we shall ommit uj dependence in E(r,o;) and x(r,i^) 
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The Fourier components of the susceptibility in the weakly deformed crystal (now depending from the 
coordinate r) can be defined according to Eq. (^) as 

Xh(r) = xf° cxp Hhu(r)] e~ w ^. (10) 
We shall look for the solution of equation (||) in the form of the expansion analogous to the Bloch waves, 

E(r)=^E h (r)exp(ik h r), (11) 

h 

k, l= k +h. (12) 

Here ko and k^ are the incident and diffracted wave vectors and the sum has to be taken over all reciprocal 
lattice vectors h. In the case of the weakly deformed crystal, when inequality (||) is satisfied the amplitudes 
E/j (r) in the expansion ( |ll|) are slowly varying functions of coordinate (on the contrary to the Bloch waves 
in a perfect crystal, when they does not depend on r). This amplitudes vary significantly on the distances 
much bigger, then the X-ray wavelengths (of the order of extinction length L ex that will be defined later). 
Therefore, if we neglect the second derivatives of E^ (r) we can obtain from (||) the following set of equations, 



ds h 
where 



? E h (r) = ^J2 \Xhh'(r) ~ a h ,5 hh ,} E fe ,(r), (13) 



h' 



kn d , > k^ , 

ah = -£T ; di-r {ShV);Sh = wy (14) 



In Eq. 



Xhh>(r) = xZl exp [-i (h - h') u(r)] e~ w ^\ (15) 

and both the displacement field u(r) and the Dcbyc-Waller factor W(r) arc slowly varying functions of 
coordinate r. 

Equations (|l^) arc the general case of the so-called Takagi-Taupin (TT) equations J|7], ^8], ^9) for the 
determination of the amplitudes of the wave fields in the weakly deformed crystals [J In the limit of a perfect 
crystal we have in Eq. ( [l5| ) for functions u(r) = and e~ w ^ = 1 and in this case Eqs. ( fl3| ) will define the 
wave field in an ideal crystal lattice. 

Taking into account that susceptibility of the crystals in x-ray range of wavelengths is small (xh ~ 
10 -5 -j- 10~ 6 ) it is possible to remain in equations (|l3|) only the waves E^(r) satisfying Bragg condition, 



For definite directions of the incident x-rays condition (|l^) can be fulfilled simultaneously for a number 
of waves. It is so-called case of multiple wave diffraction (see for e.g. book 38 and review papers |39|, p0|). 

From the other hand it is possible to find directions for which the condition (|l^) can be fulfilled only for 
one reciprocal lattice vector h, it is so-called case of the two- wave diffraction. Further we shall restrict oursclf 
only for this case. Moreover we shall consider, that the deformation field in a crystal u(z) and the static 
Debye-Waller factor e~ w ^ depend only from one coordinate z, which is the distance from the entrance 
surface to the depth of the crystal and we shall neglect its dependence along the surface. 

The x-ray amplitude of the total wave field in such a crystal in the two- wave approximation is the coherent 



superposition of the incident and diffracted waves and according to (11) is given by 



E(r) = Y, [eosE s(z)e lk ° r + e hs E hs (z)e lk ^] , (17) 



2 In the case of the strong deformation fields, when condition (g) is not satisfied, second derivatives of the amplitudes E^(r) 
also have to be taken into account pi) . 

3 In the case of the crystal with statistically distributed defects another approach of so-called statistical dynamical theory 
was elaborated (see for review Pq | and also papers |36| 
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where eo and &h are the unit polarization vectors and s is the polarization index. In the x-ray diffraction 
theory they are usually defined (see Fig. 3) respectively to the so-called scattering plane i.e. the plane 
containing the vectors ko and k/j. Polarization vectors normal to the scattering plane are called er— polarized 
(in the case of two- wave diffraction eoo-||e/ l<T ) and polarization vectors lying in the scattering plane are called 
7T— polarized (in this case polarization vectors eo^ and are misaligned by the angle 26b)- 

Now directly from the TT equations (|l3|) for the scalar amplitudes Eq(z), Eh(z) and for the fixed 
polarization s we have, 

XooEos(z) + X ohCe^- w ^E hs (z) , 

(Xhh - a) E hs (z) + XhoCe-^- w ^Eos(z)} . (18) 



dE 0s (z) 


in 


dz 


A70 


dE hs (z) 




dz 





Here f(z) = hu(z); 70.^ = cos(n • ko ; /t) are the direction cosines, n is the inward normal to the entrance 
surface of the crystal and A is the wavelength of radiation. For Bragg geometry of diffraction 70 > 0, 7/j < 
and for the Laue diffraction 70 > 0,7/! > 0. The parameter a is characterizing the deviation of the wave 
vector ko from the exact Bragg condition, 

a=^ 1 ^*-2 S m26 B (6-6 B ), (19) 
fc 

where 9b is the Bragg angle; C is the polarization factor defined as, 

q_\ 1, cr — polarization . . 

1^ cos 2#b.7t — polarization 

In most of the situations considering only the strongest elastic scattering and the photoelectric scattering 
process in dipole approximation we have for the Fourier components of the susceptibility in Eq. (p"s|): 

Xoo = Xhh = Xo, Xoh = X-h = Xh and Xho = Xh- 

Takagi-Taupin equations ( p"8] ) have to be supplemented by the boundary conditions, that for a crystal of 
thickness L have the following form for the different geometries of diffraction 

E Os (z)\ z=o = E^,E hs (z)\ z=L = (21) 

for Bragg geometry and 

E Os (z)\ z=o =E^\E hs (z)\ z=o = (22) 

for Laue geometry. 

Having in mind further applications it is convenient to transform from the set of equations ( |l8|) to a 
single nonlinear equation in the form of the Rikatti equation for the amplitude function 



where (3 = "fo/\lh\ for Bragg and (3 = 70/7/1 for Laue geometries of diffraction and Y = yXhlXK = l Y I 
exp(i^y) (for centrosymmetric crystal with monoatomic lattice | Y \= 1, $y = 0). Substituting new function 
R(z, 9) (||) into (|l|) we obtain 

T iL ex ^^ = 2[-y(9) - ty + y v ]R(z, 9) + C 1 [1± i? 2 (z, 9)}. (24) 

Here the upper sign correspond to Bragg diffraction and the lower one for the Laue. We also have introduced 
the following notations: the angular deviation from the exact Bragg position is measured by the dimcnsionlcss 
parameter, 

,(,^^MlM ± ^, (25) 
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parameters 



Xoi(l±/?) , , , ,L ex dip{z) 
W3 = ± n mv and Vv\ z ) = ± - 



2^X r "W"-/-- 2 (26) 

define attenuation of x-rays due to the photoelectric absorption and the shift of the Bragg position due to 
deformation in a crystal; 



C 1 = C{l-i P )e- w ^\p = -^-- 



Xi 
X 



L ex is an extinction length defined as 



A70 



(27) 



(28) 



Here we have also introduced the following parameters X r = Re^/XhXh and X = Im^XhXh- Now boundary 
conditions for equation (|24|) are defined on one surface. For the Bragg case of diffraction we have R(z) \ z =l = 
and R(z)\z=o = for Laue case. 

The reflectivity is usually defined for Bragg case as 



P R {6) = (l/0)\E h (O,6)/E o (O,i 



(29) 



now has the following form. 



P R {0) = \Y-R(p,0)\' 



(30) 



It is easy to obtain solutions of the equation ( p4[ ) in the case of a perfect thick crystal (iM)L » 70, where 
Ho is a normal absorption coefficient defined as /io = kxoi)- In this case <p(z) — 0,e _M/ ^' = 1 and Eq. ( |24| ) 
reduces to an equation with constant coefficients. So, for thick perfect crystal solution does not depend on 
the thickness of a crystal, that is we have dR/dz — 0. Now from (|i|) for Bragg case we obtain directly 



Ro(0) 



1 

"Ci 



{-y - iyo) + \J(y + m) 2 -cf 



(31) 



where for the square root it is chosen the branch with the positive imaginary part. 

For the amplitude of the refracted wave E 0s (z,9) we have directly from the TT equations ( |l8| ) (and 
taking into account definition (p3|)) 



dE 0s (z,6) 
dz 



™XQ ■ Cl r>, a , 

t—R(z,9) 

A7o Lex 



E 0s (z,9). 



Formal solution of this equation can be written in the following form 

z 

inXo ■ 1 



E 0s (z,9)=E^ 



exp 



z-i / dz'CxRiz'^ 

A70 L e 



and we have for the intensity of the incident wave, 



(32) 



(33) 



I (z,9) = \E 0s (z,6)\ 2 =lt ) exp 



Mo . 2 
— z H 1m 

7o L ex 



z 

J dz'C x R{z', 



In the case of a perfect crystal, R(z,9) = Ro(9) and we have from (|34|), 



I (z,6)=4 in) exp 



7o 



= 4 in) exp l-^-z + ?Llm [CM9)} 

I 70 ^ex 

'We want to note, that our choice of extinction length differ from commonly used by the factor it. 



(34) 



(35) 
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where fii n (9) is an interference absorption coefficient. This expression takes not only into account normal 
attenuation of x-rays out of the angular region of the dynamical diffraction (y >> 1) 

I (z,6)=lt ) exp{-^z\ (36) 



To 

but also takes into account a dynamical " extinction" effect coming from the multiple scattering of x-rays on 
atomic planes in the narrow angular region of the dynamical diffraction |[], 0, || . In the region of the total 
reflection for y ~ 0, we obtain from ( p5| ) 

/oMW^cxpj-gz}. (37) 

Here we have taken into account also that yo << 1 and [Iqz << z/L ex . From this expression we can see that 
for the angular position y = x-rays are effectively attenuated on the typical distances z ~ L ex that for the 
energies E ~ 1 lOkeV are of the order of microns and are much smaller then normal attenuation distances 
z ~ 70/^0 that for the same energies can be of the order of tenth and hundreds of microns (see e.g. 0]). 

As we can see from the expression ( |37| ) extinction depth L ex is one of the important parameters of the 
theory that give an effective attenuation distance for x-rays while the dynamical diffraction. In our further 
treatment all other distances will be compared with L ex . 

Here we want to make several remarks. The amplitudes Eq s {z, 9) and Ehs(z, 9) in TT equations ( |l8| ) are 
complex numbers with its amplitude and phase. Due to the fact that the dynamical scattering is a coherent 
scattering process this two amplitudes are connected with each other and, for example, in the case of a 
perfect crystal on its surface we have from ( fj3| ) for the ratio of these amplitudes on the surface of the crystal 



E hs (z,6) 



E 0s {z,( 



z=0 ROWI 



where Rq is defined in (|3l|). 

Typical behavior of the reflectivity Pr(9) and of the phase a(9) in the diffraction region is shown on Fig. 
4. In this small angular region typically of several arcsec the reflectivity Pr{9) is of the order of unity and 
the phase a(9) of the wave field changes from — n to 0. Just this fast change of the phase makes x-ray 
standing wave method so sensitive to any additional phase shifts. 

2.2 Susceptibilities 

The Fourier components of the susceptibility xo and Xh (see expansion (Jfy) are in general complex valued 

Xh = Xhr + iXhi- (39) 

The real part Xhr correspond to elastic scattering of x-rays and imaginary part Xhi accounts for absorption 
effects. The values of Xhr and Xhi are calculated from quantum mechanics (see Fig. 5, where the values of 
Xor and xoi are calculated for Si and Ge for different energies) and for crystals without center of symmetry 
may themselves be complex ||. For hard x-ray energy range (E ~ 1 lOkeV) Xhr is negative and for the 
most of elements is of the order of 10~ 6 . It is convenient to present it in the following form 0], 



Xhr 



F hr , F hr = J2 (fj + A/i) e~ Wi e~ iilpi ■ (40) 
j 

Here rg = e 2 /mc 2 = 2.818 • 10~ 15 m is the classical electron radius, il is the unit cell volume and Fh r is 
the structure factor for the reciprocal lattice vector h. Expression (^) is written for an arbitrary unit cell 
of a crystal, summation is made over all atoms of the unit cell, Pj is the coordinate of the j— th atom in 

5 Notc, that wc have defined the E-field as (a) E 0ih e l ^o,h^-^) ( scc Eq. (2.15)), whereas frequently (b) _E 0ih e i(a "~ k o.'' r) is 
used. However, this only introduces different phase convention a a = —a^ if we denote the phase a resulting from the case (a) 
and (b) with a a and a b . 
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a unit cell; e w j is the thermal Debye- Waller factor that takes into account the attenuation of the elastic 
scattering of x-rays due to a thermal vibrations of the atoms. In equation ( [Ili| ) 

/ 3 -(h) = J n(r)e^dr (41) 

is an atomic scattering factor for the j— th atom in a unit cell. It is determined by the electron density n(r) 
in an atom and Afj is an account for the dispersion corrections to an atomic scattering factor. The values 
of this parameters are tabulated in International Tables for X-ray Crystallography H . 

As it was already mentioned above, the imaginary part of the susceptibility xm takes into account 
absorption effects. For hard x-rays (E ~ 1 -j- lOkeV) its value (see Fig. 5) is two orders of magnitude smaller 
then the real part \hr (Xhi ~ 10 -7 -j- 10~ 8 ) and, for our choice of the phase in the plane wave (|TT|) , it is 
positive. It can be shown p3| , [42| that in general case the imaginary part of susceptibility Xhi contains 
contributions from all the inelastic processes: the photoelectric absorption, Compton scattering and thermal 
diffuse scattering 

Xhi = Xhi(Ph) + Xhi(CS) + Xhi{TDS). (42) 

The imaginary part of the Fourier component of the susceptibility Xhi m a crystal in analogy to (^) can 
also be presented as a sum of contributions of different atoms 



j 



where <jj are the cross sections of the different inelastic processes for the j— th atom in the unit cell and their 
values can be obtained from j|, f43[ . 

As it was pointed out previously in general case the susceptibility of a crystal is a tensor and has a 
non-local character. Being interested in diffraction and taking into account relationships between the values 
of the cross sections of the different processes (sec Fig. 6) 

OT » OPh » 0"C > &TDS, (44) 

where ctt, uph, &c, &tds are the cross sections of the elastic Thompson scattering, photoelectric absorption, 
Compton scattering and thermal diffuse scattering we can neglect in (|l8| ) small non-local corrections to Xh 
and treat susceptibilities as scalar values. However analysing the yield of the secondary radiation while the 
dynamical diffraction of x-rays this corrections may be essential and can not be neglected. For example, while 
considering Compton and thermal diffuse scattering it is necessary to account tensor character of Xhi and the 



angular dependence of the corresponding cross sections of the inelastic scattering (see for e.g. |4J, |45j, |46|). 
Different situation is realized for practically valuable case of fluorescence radiation and photoelectron emission 
(according to ( piij) it is the main inelastic process). So far as these processes are caused by photoelectron 
absorption, the total value of each in dipole approximation does not depend from the direction of propagation 
of the radiation in an isotropic crystal the imaginary part of the susceptibilities can be treated as scalar values 
and without angular dependence. Further, if not mentioned specially, we shall consider mainly this case. 

Takagi-Taupin equations ( |l8| ) were obtained in the dipole approximation. Small quadrupole corrections 
in imaginary part of Xhi, if necessary, can be also taken into account. They will bring to renormalization of 
the polarization factor C, that can become essential for scattering near adsorption edges and backscattering 
(see for details @ pi |§). 
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3 Theory of x-ray standing waves in a real crystal (general ap- 
proach) 

In this Chapter we shall obtain, using the approach of Afanasev and Kohn p^ |, the general expression for 
the yield of the secondary radiation in the case of the dynamical diffraction of X-rays from the deformed 
crystal lattice. The amplitudes of the waves Eq(z,9) and Eh(z,9) in such a crystal can be obtained from 
the TT equations (|TJ) (we shall consider deformations that depend only from z). 

To find the intensity of the secondary radiation yield at the depth z in a crystal, one must determine 
the number of absorbed quanta in a layer with the thickness dz per unit area and unit time (Fig. 7). It is 
proportional to the loss of the energy field in this layer. From the equation of the field energy balance we 
have for the number of absorbed quanta 



dN(z) 
dz 



-divS(z) 



dz 



(45) 



where S(z) is the energy flow (Poynting vector), averaged over the time period of the field oscillations and 
over the elementary cell of the crystal. In (|4^) we have taken into account that S(z) depend only from z. 
According to the definition of the energy flow vector 



S(z) 



8tt 



[s Q \E Q {z)\ 2 + s h \E h (z)\' 



(46) 



where So = ko/|fco| and = k^/|fch| are the unit x-ray propagation vectors, we can obtain for the number 
of absorbed quanta, 



dN(z,6) 



dz %-kTluj 
Taking into account TT equations (E^) we have 



7o- 



d\E (z,i 



dz 



lh - 



d\E h {z,i 
dz 



(47) 



dN m (z,0) 
dz 



{E*(z,6) X ooA m ) E o(z,Q) 



ck 

EZ(z,9) X hh,i(m)E h (z,0) + 
2Re{El(z, 6)xohA m ) E h(z, 9) exp(iip(z) 



(48) 



W(z))}}. 



The number of absorbed quanta is determined only by the imaginary part of the susceptibility Xhi that 
is account for absorption effects. According to (42) it is possible to separate the influence of the different 
processes into the yield of the secondary radiation. The index m introduced in (fi"8|), characterizes this 
contribution of a certain secondary process, which is under investigation. 

The total number of the secondary quanta emitted from the crystal is equal to 



N m (9) = / dzP™(z) 



d N m (z,6) 
dz 



(49) 



where P™(z) is the probability function of the yield of the secondary radiation of the type m from the depth 



Equations @[49|) are general and give the solution for the problem of the angular dependence of the 
secondary radiation yield, when x-ray standing wave exist in a crystal. They are valid for any type of inelastic 
process such as photoeffect, fluorescence radiation, Compton scattering and thermal diffuse scattering. They 
can be applied as well for investigation of secondary electrons, i.e. Auger electrons and electrons ejected 
due to absorption of fluorescence radiation. One must only define the values of the Fourier components 
of the susceptibility Xoo,i>Xhh,ii an d Xoh,i appropriately, as well as the probability function P™(z). The 
amplitudes E$(z,9) and Eh(z,9), naturally, does not depend on the type of the inelastic process that is 
experimentally registered, but are determined only by the diffraction process on a real crystal. If, for 
example, the deformation field u(z) and the level of amorphization W(z) are known, then the amplitudes 
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E a (z, 9) and Eh{z, 6) can be obtained directly from the TT equations (|l8|). On Fig. 8 results of calculations 
of the reflectivity curves and photoeffect yield for the silicon crystal with the known profile of deformation 
(also shown on Fig. 8) are presented. The series of curves correspond to both the entire layer (upper curves) 
and to parts of it. 

In most of the applications of XSW method it is assumed, that the yield of the inelastic process under 
investigation is completely determined by the intensity of the wavcficld at the atoms positions. However, 
according to result summarized in Eqs. (fl8|-|49|), it is not always fulfilled. In fact there are two main effects 
that are taken into account. This is first of all the deformation of a crystal lattice described by an additional 
phase factor <p(z) (due to the displacement of atomic planes) and the static Debye- Waller factor e~ w ^ (due 
to the random displacements of atoms) in the third term of ( |48| ) . The second effect is coming from non-locality 
of some of the inelastic processes. This can be, for example, effects of higher order multipole interactions for 
the photoeffect processes [S, Q or non-local character of such processes as Compton scattering or thermal 
diffuse scattering @ f£j| [Sf. 

Being interested in future mainly by the fluorescence and photoelectron yield in dipole approximation we 
can write Eq. (|i^ ) in the following form (further we will omit index m and assume, that xoo,i = Xhh.i = Xoi 
and xoh,i =Xhi), 

dN iz ,e) _ c/xo |SoM)|2{1 , , 



dz fyrhu' \E h (z,6)\ 2 

+^^§^4^ W{z) }h (50) 



where 



2>. 



wT-ih Pi 



e j e 



e Ti = ^ = ^ . (51) 

XOi E 3 

3 

Now substituting into (|50| ) the expression for the amplitude R(z, 6*)(^3|) we have for the normalized intensity 
yield of the secondary process 

oo 

<0) = ^ = J—fdzP yi (z)\E (z,6)\ 2 {l + l3\Y\ 2 \R(z,e)\ 2 + 
/(oo) /(oo) J 



+2^e- w ^CRe[e 1I YR(z,e)}}, (52) 
where intensities are usually normalized by their values far from the region of the Bragg diffraction /(oo) = 

oo 

fdzP™(z)\E (z,oo)\\ 

Equation ( |52"| ) together with equations (^4|, ^) and (|35| ) completely determine the scheme of calculation 
of the angular dependence of the yield of the secondary radiation in the most general case under the condition 
that a plane wave is incident on the crystal. In a real experimental situation the experiment is performed in 
a double-crystal scheme with a first crystal-monochromator. In this case for the comparison of theoretical 
calculations with the experimental results the convolution between the curve of the secondary radiation yield 
and the reflectivity curve of the monochromator crystal has to be calculated. If an asymmetric reflection 
is used in both crystals, i.e. asymmetry factor (3 in the sample crystal and the asymmetry factor of the 
monochromator crystal fli do not equal to unity, then we have for the convoluted curve of the SR yield 



K(6) = 4P- = . (53) 

s — oo 

Here 8 is the angle between the reflecting planes of the crystal-monochromator and the sample crystal and 
summation in Eq. (J53|) is performed over the different polarization states. Convolution with the reflectivity 



curve (30) of the sample Pr(A6) is defined in the same way. 
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Finally, in this Section we have formulated the main equations for the SR yield excited by XSW while 
the dynamical diffraction of x-rays in real crystals. In the remainder of the work we will analyse different 
physical applications to this general formalism. 
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4 XSW in a perfect crystal 



In this Chapter wc shall consider an effects of registration of the different inelastic processes in a simple case 
of a perfect crystal. Though it is the most simple case of the real crystal, the main peculiarities of the XSW 
field and SR yield can be already revealed and understood in this case. In the end of the Chapter the case 
of amorphous layer on the top of the perfect one is shortly discussed as well. 

4.1 XSW with a big and small depth of yield (extinction effect) 

In the case of the perfect crystal the amplitude R(z, 9) = Rq(9) (see Eq. (|3lj)) and it does not depend from 
the coordinate z. From (52) we obtain for the wave held intensity, 



1(9) = {1 + (3\Y\ 2 \R (9)\ 2 + 2^CRe[e JI YR (9)}} / dzP yi (z)\E (z, 6)\ 2 . (54) 



In the limit when we can neglect the angular dependence in the last integral and in addition if approx- 
imation e-f- = Xj[i/xoi = 1 is valid we are coming to the well known expression for the wave field intensity 
(|J) . From the above expression we can see that the shape of the curve of the SR from the perfect crystal 
is mainly determined by two factors. First of all it depend on the depth of yield L yi of the SR that is 
determined by the probability yield function P y i(z) and, secondly, in the case of detecting SR process from 
the multicomponent crystal its shape essentially depend from the complex factor s-r- In the beginning we 
shall consider effects of the depth of yield. 

For the future analysis it is convenient to take the probability yield function P V i(z) in the form of an 
exponential function: 

P yi (z) = exp(- M ^). (55) 

This form of the probability function is exact for the yield of the fluorescence radiation with p e JJ = p y l/jfi, 

where \i yi is an attenuation coefficient of the specific fluorescence line, that is measured in experiment and 7/; 
is the cosine of the exit angle of this fluorescence radiation. For the case of photoeffect integrated over all the 
directions of the photoelectron yield this form of P y i(z) with (j^y = 2.3/L^, where L v ^ is an average escape 
depth of electrons is a good approximation |5l| to probability yield function obtained from the Monte-Carlo 



simulations 1 52 , B3| . For the escape depth of the electrons with an initial energy Ei (in keV ) the following 



approximation formula can be used 

l$(Ei) = 780E1/PHE.JEO), (56) 

where p is the density of the material in g ■ cm~ 3 and Eq w 0.39keV. 

Taking now into account that according to Eq. (35) for the perfect crystal \Eq(z, 9)\ 2 = \Eq' | 2 exp(— /ii„ (9)z/ r yo), 



we obtain from (154) for the yield of the secondary radiation from a perfect crystal 



1(9) = J- ^ , (57) 

4 f + ^(o)ho 



where 



I sw (9) = 1 + I3\Y\ 2 \R (9)\ 2 + 2^3CRe[e^YR (9)}. (58) 

Now we can easily analyse different limits of the depth of yield parameter on the angular dependence 
of the SR curve. For example, in the limit when the escape depth of the secondary radiation L y i is much 
smaller then the minimum penetration depth L ex of the standing wave field into a crystal (L y i << L ex ) 
p e JJ » p,in(9) and from Eq. (|57j), we obtain 



1 p m (9) 
7o p e y { f 



(59) 
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The angular dependence of the SR yield in this limit is mainly determined by the intensity variation 



of the standing wave held through the atomic planes (Eq.(58)). Its shape for the monoatomic crystal will 
totally coincide with the intensity variation of the wave field (|). At the same time, following Eq. ( |59| ) the 
maximum change in the shape of the intensity curve is due to an extinction effect that can be seen only in 
the central region of the total reflection. It leads to a weak variation in the slope of the linear part, i.e., to a 
weak decline in the intensity yield at this angles. It is interesting to note here, that this actual form of the 
standing wave curve contain information about the real escape depth of the secondary radiation, that can 
be obtained by fitting experimental data to theoretical calculations in the form of equation ( |59| ) . This effects 
were observed experimentally in the case of the fluorescence radiation |54|] (change of the depth of yield L y i 
was obtained by the change of the exit angle of the fluorescence yield) and in the case of photoeffect |5q ] 
(see Fig. 9). 

Now we shall consider opposite limit of the big depth of yield (L y i » L ex ). This is the typical situation 
in the case of measuring fluorescence radiation from the atoms of a crystal lattice or for the measuring of 
the inelastic scattering on thermal phonons. For the analysis of this limit it is useful to use in Eq. (p4}) the 
following expression (see for e.g.Q), 

Mo {l + 0Pn(e) + 2^/]3Re [ejR (9)] } = ^(6) [1 - P R (6)} (60) 

and now we obtain for the yield of the SR, 



1(9) = [1 - P R (9)] e/ f" (g)/ f° - [1 - Pr(0)] 



eff 



(61) 



In the case of the fluorescence radiation we have for the intensity 

P l {6) = [1 - P R (9)] [l - (7o/7/OO^/)M0)) 

As we can see from Eq. ( |6l"| ) now in the limit of the big depths of yield /j, y i << Hi n (&) the shape of the 
SR curve has the form of the reverse curve of the reflectivity 1(9) ~ [1 — P R (9)]. This result has a simple 
physical explanation. Really, in the angular region of the dynamical scattering x-rays does not penetrate 
deeper than the extinction depth L ex , so, SR can be excited only from this depths. Due to the law of the 
energy conservation the yield of the SR has to be equal to Eq. (|6l|). However, if the term lo^JJ / ^in(9) in 
expansion ( pl| ) is becoming comparable with unity, then on the curves of the secondary radiation one can 
see small asymmetry due to the behavior of /ij„ (9) . This behavior of the fluorescence radiation was for the 
first time measured and understood in the pioneer works of B. Batterman [?J and then repeated in many 
other works (see for e.g. Fig. 10). 

It is clear, that using the depth of yield L y i as a parameter we will obtain a number of curves that lie 
between two limiting cases described by Eq. ( p9| ) and (|6l|). On Fig. 11 the calculated curves of the angular 
dependence of the SR yield from a perfect Si crystal are presented. For calculations it was used the case of 
(400) diffraction of CuK a radiation (L ex = 3.60/im) for different values of the parameter L y j. 



4.2 Multicomponent crystals 

In a perfect crystals with the complicated elementary cell containing different type of atoms XSW method give 
a unique possibility to investigate position and degree of disorder of different type of atoms. Experimentally 
the most effective way to do it is to register the characteristic fluorescence radiation from different atoms. 
The nodes and antinodes of x-ray standing wave are located in the different way for different sublattices of 
the crystal (on Fig. 12a (111) diffraction planes of GaAs crystal are shown). In this situation the angular 
dependence of the fluorescence yield for different type of atoms will have the different shape and also will 
differ from the typical curves of the monoatomic crystals (see for e.g. Fig. 12b for the same case of GaAs 
crystal) . 

Peculiarities of the angular dependence of the fluorescence yield in a multicomponent perfect crystal 
are determined in fact solely by the factor e-r in (Q). If we are interested by the yield of the fluorescence 



16 



radiation from the atom of the sort a from a multicomponent crystal then we have for the factor e-^ in ( |54| ) 
(see also @), 

, _ N E^- ff '<%«(h) 

where c" is the concentration of the atoms of the sort a in the sublattice j, a a (O)and <r (h)are the cross sec- 
tions of the corresponding processes, e~ w i are the Debye- Waller factors (here they are the sum of thermal 
and static displacements) and 5?(h) = J] e 1 p 3' are the structure factors corresponding to the positions of the 



3 



atoms of the sort a in the unit cell. The total cross sections of the photoexcitation in dipole approximation 
are isotropic and we have for cr a (h)/er a (0) = 1, however if quadrupole contributions are valuable we have for 
the ratio of cross sections in Eq. < p2jj 

(63) 



o-„(0) a a V C 

where a a = c<? + ^a^a an d o® are the total, dipole and quadrupole cross sections correspondingly, C 



is the polarization coefficient (20) and parameter is equal to cos 26* s for cr— polarization and cos 46* # 
for 7T— polarization (see for details j|8l ^9|). This can be important while measuring fluorescence radiation 
near absorption edges or in a backscattering geometry, where quadrupole contribution can be essential. 
Expression ( |62| ) is simplified for the case of two-component crystals when different type of atoms occupy 
different sublatticcs (this is the case of GaAs, InSb and etc. crystals). In this case we have from Eq. ( |62] ) 

e|~ e- w "^S a (h) = e - w "W | s a (h) | e lv,a(h) , (64) 

where dipole approximation is assumed. Structure factors Sj(h) are complex quantities with their amplitude 
| S' a (h) | and phase f a (h). Substituting this values of Sria) into expression ( |5S| ) we obtain for the intensity 
of standing wave on the positions of the atoms of the sort a, 

I sw (9) = 1 + f3\Y\ 2 \R Q (9)\ 2 + \ YR Q {9) \ F?(h) cos [a{6) + ^ a (h) + $ r ] , (65) 

where 

F?(h) = e - wa ^ | S a th) | . (66) 

This is an important result. According to Eq.(^5^the angular dependence of the SR yield directly depend 
on the phase ip a (h.) of the structure factor S^h). Moreover, due to the fact that this phase enter the 
interference term in ( |6"5| ) the angular dependence of the field intensity I sw (6) is very sensitive to the value 
of this phase. So, measuring this angular dependence one can determine with high accuracy the phase 
and the amplitude | S* a (h) | of the structure factor for different sublattices (or different sort of atoms) in 
multicomponent crystals. In this way the h— th Fourier component of the structure factor can be totally 
determined J58|. This is illustrated for the case of the GaAs crystal on Fig. 12. Of course according to 
our previous discussion this initial curve will be modified when the final depth of yield of the fluorescence 
radiation (Eq. ([37])) will be taken into account. 

This idea of measuring the SR (fluorescence or photoelectron yield) in the XSW field in multicomponent 
perfect crystals was successfully realized in a number of experiments. For example, the polarity of GaP 
crystals was obtained in experiments ]59| while monitoring fluorescence radiation and in experiments ]6C| | 
GaL photoclcctrons were measured using a cylindrical energy analyzer in a high- vacuum chamber. Due to a 
high sensitivity of XSW method it has become possible to measure the change of the phase of the structure 



factor as a function of energy near absorption edges in noncentrosymmetric crystals (see for e.g. 57, |6l| ] 
and also theoretical paper (6^]). In the paper |Q it was demonstrated the possibility of determination of 
the positions of the Ga and Gd atoms in the unit cell (it contain 160 atoms) of the perfect garnet crystals 
while monitoring characteristic fluorescence radiation in the XSW field. In a recent papers the fluorescence 
radiation from different atoms in the unit cell of High-T c crystal was measured |34j, positions of Cd, Zn, Se 
and Te atoms in Cdi- x Zn x Se y Tei- y single crystals were obtained ]6^ j. 
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Analysis of the photoelectron yield (especially if Auger electrons are monitored) while scanning XSW in 
multicomponent crystals is more complicated comparing to the fluorescence yield. Really, even if detectors 
with high energy resolution are available then the contribution to the total photoelectron yield in the reg- 
istered energy range is in general the sum of contributions of primary photoclcctrons ejected from different 



atomic subshclls and from different sort of atoms. Finally we have 12 



Y el {E, 6) = J2 9NaY Na {E,6), (67) 

N,a 

where for the case of perfect crystals Yj^ a (E, 9) is determined by Eq. ( |54| ) with the factor e-^ = ef-(TV) equal 
to that of ( p2| ) and with P y i(z) = Pjsj a (z, E), that determines now the probability of the photoelectron escape 
with energy E ejected at the depth z from the TV— th subshell of the a— th sort of atom. Parameter g^a 
determines the fraction of electrons ejected from the TV— th subshell of the a— th sort of atom to the total 
number of such electrons, 

n a a Na (0) 

9Na - j^r, (08) 

N,a 

where n a is a number of atoms of the a— th sort in a unit cell. Taking into account, that typically the 
photoelectron yield depth of yield is much smaller, then the extinction depth (L yi << L ex ) we can use for 
the photoelecton yield expression ( |58| ) with the parameter e-g equal to, 



Y,9Nael(N)P Na {E), (69) 



JV.a 



where P Na (E) = J P Na (z,E)dz. 
o 

4.3 Crystals with an amorphous surface layer 

We shall assume now, that the top of the crystal contain thin amorphous layer with the thickness L am . We 
shall also assume, that the depth of yield of the SR is much smaller, then the extinction depth, but bigger 
then the thickness of amorphous layer L am < L V i << L ex (this can be the photoelectron yield for example). 
It is clear, that in amorphous layer there is no interference between incident and diffracted beam and we can 
neglect small attenuation of x-rays in this thin layer. Taking all this into account and performing integration 



in (52) separately in amorphous and in the perfect part of a crystal we obtain, 

«(*) = P am n am (9) + (1 - Pam )K* d (9), (70) 

where 

K am {9) = l + I3P R {9) 1 (71) 



K ld (9) = I sw {9) = 1 + I3P R {9) + 2^pCRe[e~YR {6)\ (72) 
and p am is amorphous fraction of the crystal defined as 



Pam = J Pyi(z)dz. (73) 


oo 

In (|7C|-f73"|) we are assuming, that the depth of yield function is normalized to unity J P y i(z)dz = 1. 



From expressions (|70|-[72[) we see, that in general for different thickness of amorphous layers we obtain 
different shape of the curves of the photoelectron yield all lying between two curves K ld (9) for perfect crystal 



18 



and K am (9) for amorphous layer L am > L y i. Such measurements were made in |6(], |67J and are presented on 
Fig. 13. 

In general escape depth of electrons L V i(Ei) from the crystal depend from the initial energy Ei (see for 
e.g. Eq. (56)). Assuming now exponential probability yield function in the form (pq) 



1 

Lyi(Ei) Lyi(Ei) 



we obtain for the amorphous fraction of the crystal p am 

Pam = 1 - exp[-L am /Lyi(Ei)]. (74) 
From this expression we have for the ratio L am / L y i(Ei) 

L am /L yi (Ei) = -ln(l -Pam)- (75) 

Now, if one of the parameters: thickness of amorphous layer L am , or escape depth of electrons L V i(Ei) is 
known, then another one can be obtained from Eq. J7q). 

This idea was realized in the paper |5J| (see Fig. 14), when the photoelectron yield for different group of 
electrons with different loss of energy was measured with the low-resolution gas-proportional counter. In this 
experiments escape depth L yi = 1.2 fim for the initial energy of the electrons Ej = 13.3 keV in Si crystal 
with amorphous Si0 2 layer was obtained. 
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5 Bicrystal model (Bragg geometry) 



5.1 Theory 

The calculation scheme of the angular dependence of the SR yield k(6) for the Bragg geometry of x-ray 
diffraction in the crystals with the deformed surface layer was developed in J50]| . In general case when the 
deformation profile of the crystal lattice has an arbitrary profile the main problem is to solve nonlinear 
differential equation (p3) for the amplitude R{z,9) (p3|). For the arbitrary dependence of the functions 
tp(z) = hu(z) and e _v ^' from the depth z Eq. (|24|) and integrals in ( |34] ) and (|5|) can be calculated only 
numerically (see for e.g. Fig. 8). 

In this Chapter we will discuss the simplest model of the deformed crystal, so-called bicrystal model, 
that allows to find analytical solution. In the frame of this model crystal consists of two parts a thick 
perfect substrate and the deformed layer of thickness Ld with the linear dependence of the deformation field 
u(z) = (Ad/d)z, where parameter Ad/d is a constant difference of the interplanar distance in the layer 
comparing to that in the substrate (Fig. 15). In this layer parameter y v {z) = y v = n(L ex /d)(Ad/d) ( p6| ) 
and static Debye-Waller factor e~ w W — e~ w ° has constant values. So, in this model deformed layer has 
a different interplancr distance comparing to that of the perfect substrate. It is additionally partially and 
uniformly amorphysized and sharp transition (in fact step function) between the layer and the substrate is 
assumed. Though this model is very simplified it turned out to be very practical in the analysis of the SR 
yield in a number of experiments. 

Analytical solution for the angular dependence of the SR yield k(9) and reflectivity Pr{9) for Bragg 
diffraction and bicrystal model was obtained for the first time in |5l| . This result is a particular case of a 
more general approach of a multilayer crystal consisting of a number of layers with the different parameters 
that will be discussed in details in Chapter VI I (se e also Appendix). According to this approach SR yield 
for bicrystal model can be obtained from Eq. (|l32| ) with the number of the layers N = 2 

<oM = K a Ld {°) + \T(L d ,6)\ 2 e-^ L «K a (6). (76) 
Here, normalized SR yield nf ot (9) is a sum of two contributions: from the deformed layer K a h (9) and 



a perfect substrate Kq(9). This functions are defined in Eq.(129) and transmission coefficient |T(L<j, /1> 



can be determined from (122h. Deformed layer is characterized by constant values of deformation Ad/d, 
static Debye-Waller factor e~ w ° , coherent fraction F£ and coherent position P". For the perfect substrate 
Ad/d = 0, e~ w ° = 1 and, in general, it can has the values of coherent fraction F" and coherent position P c a 
different from that of the layer. The angular dependence of the reflectivity Pr{6) for Bragg diffraction is 



obtained from (|30|) and (114-116) with the amplitude Rd„ = Ro(9), where Rq(9) is the value of the amplitude 
for the perfect substrate defined in (|3l|). 

We want to note here that obtained result is exact and is valid for any relationship between the values 
of Ld, L ex and L y i. It is also valid for the analysis of the angular dependence of the photoclcctron yield 
or fluorescence radiation in the XSW field. Calculations performed with an exact expression ( |76| ) give 
possibility to take into account minor effects of extinction and variation of phase on the depth of yield of the 
SR. The convenience of the analytical approach in comparison with the direct numerical calculations lies in 
the possibility of explicitly separating the dependence on the various parameters. 

Essential simplification of Eq. ( |76| ) can be obtained in the special case, when the following condition 
Lyi « Ld « L ex is satisfied. It can be easily fulfilled, for example, for the photoelectron yield. In this 
case in the angular region of the dynamical diffraction from the substrate \y(0)\ < 1 the SR yield (Q) can 
be presented in the following way 

k(6) « 1 + pP R d {6) + 2^f3C \R Q (9)\F?(h) cos [a Y (9) + y a h + cp ] , (77) 

where P R d {9) = \YR a (9)\ 2 is the reflectivity from a perfect substrate (|§Q; F?(h) = \Ye^\e^ w ° is the 
coherent fraction; ay{9) — a(9) + arg(y) where a(9) is the argument of the complex amplitude Rq(9); 
(fh = 2ttP^ = arg(ei-), where P° is the coherent position; and the phase ipo = y(0) = 2y ip L d / L ex gives the 
total phase shift due to deformation in the layer. In fact this result is a special case of a more general result 
obtained for the first time in fl32| and is valid for any type of deformation if the condition L y i << Ld << L ex 
is satisfied. 
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5.2 Experiment 

Now we will show on a number of examples how this theoretical approach can be applied for the analysis of 
different experiments with the use of XSW for the investigation of the real surface structure. 

For the first time it was applied for the analysis of homoepitaxial films on Si surface (^ij (see also [[38| ) . 
In this experiments photoelectron yield was measured from a set of specially prepared Si single crystals with 
the grown homoepitaxial Si films doped with Ge. The concentration of Ge varies from sample to sample 
from 3.7 x 10 19 to 1.5 x 10 20 atoms cm -3 . The film thickness was Ld = 1.5/im. The CuK a radiation and a 
(444) non-dispersive double crystal diffraction arrangement with an asymmetric monochromator were used. 
In this case L y i = 0.45^m, L ex = 10.5/im, and consequently the condition L y i « Ld « L ex is fulfilled. 
The experimental results are presented in Fig. 16. The doping with the Ge atoms changes the plane spacing 
Ad uniformly in the disturbed layer, which leads to a change of the surface displacement. Hence the value 
of the phase on the surface (po = |h|u(0) = 2ir(Ad/d)(Ld/d) also changes. This factor along with the change 
of Debye- Waller factor e~ w ° in the surface layer is responsible (according to Eq. ([77])) for the variation of 
the shape of the photoelectron yield curve. Table 1 presents the values of <po and Wq obtained by fitting of 
the experimental curves on Fig. 16 to the theoretical ones calculated using Eq. (f77|). In this experiment it 
was demonstrated for the first time high sensitivity of the XSW method for measuring small displacements 
in the surface layer (in fact, relative difference of the lattice parameter in the layer comparing to that in the 
substrate). We want to note here, that though lattice parameter changes were comparably not so big from 
sample to sample the phase difference for high order (444) reflection was of the order of ir that brought to a 
significant change of the shape of the curves from sample to sample (see Fig. 16). 

In the paper |6l| the films of Ino.5Gao.5P of different thickness (~ 0.01/im and ~ 0.6/im grown on the 
surface of GaAs (111) single crystal were investigated. The fluorescence radiation from the In and P atoms 
excited by the XSW field in the substrate (in the case of thin film) and in the film (for thick film) was 
measured (see Figs. 17, 18). 

In the case of thin film condition L vl = Ld « L ex is fulfilled and we can use Eq. (|77|) to obtain the 
phase shift ipo of the surface layer. Taking also into account relationships (^52- 64) the angular dependence of 



the fluorescence yield from the thin film can be analysed from the simple expression, 

* a iM = 1 + ijrfjp + 2C ]1^™ cos [a{6) + ^ + ' (78) 

were parameter F^(h) (coherent fraction) is equal to F®(h) = \S^\ exp(— W a (h))Fo, S% = \S%\ exp(zy^) is a 
structure factor of the atoms of the sort a, a(9) is the phase of the complex ratio Eh/ Eq (see Eq. ( p8| ) and 
Fo = e~ w ° is the average amorphization of the surface layer. 



From the fitting of the experimental curve to the theoretical one (Eq. (78)) it was obtained the value 
of Fo = 0.25 and of the phase shift tpo = 0.55 that correspond to the average deformation in the film 
Ad/d « 6 • 10~ 3 . The same approach was used in Jfof for the analysis of the film In x Ga\- x As y P \- y grown 
on the substrate InP (100). Small deformation Ad/d ~ 2.3- 10~ 4 in a thicker film Ld ~ 0.2/ito was measured 
with XSW method. 

In the case of the thick film x-ray standing wave field is formed in the film itself. Fluorescence yield 
measured in the angular position of the film maximum has different behavior for In and P atoms due to 
different position of the atoms in the unit cell. In this case the simplified approach of Eq. ( |78| ) is not valid any 
more and general approach based on Eq. ( [76] ) was used for fitting (see Fig. 18). As a result thickness of the 
film Ld = 0.62/im, deformation Ad/d w 2 • 10~ 3 and amorphization factor Fo = 0.8 were obtained. However 
as it is seen from Fig. 18(b) coincidence between the experimental and theoretical curves is not perfect, that 
can be due to simplified bicrystal model, that does not take into account transition layer between the film 
and the substrate. In this situation the theoretical model with transition layer proposed in [fziJl can be useful 
for the analysis of the angular dependence of the fluorescence yield from the geterostructures. 

Thin epitaxial films of CaF2 grown on Si (111) surface where characterized by impurity luminescence 
probes, x-ray diffractometry and x-ray standing wave technique ]72| , [7^] . Molecular beam epitaxy was used 
to grow lOnm thick films. The GaF^/Si interface was formed at 770°C. Due to different growth conditions 
strain field in each GaF^ film is different. This is well seen from the reflectivity (222) curves measured from 
the different samples (Fig. 19(a)). The angular dependence of the CaK a fluorescence excited by the XSW 
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field in the Si (111) substrate was also measured (Fig. 19(b)). In this case we again have the situation, 
when L y i = Ld « L ex , and expression ( |78| ) can be applied for the analysis. Finally the strain field in the 
film was obtained independently from the photolumincscence study and x-ray rocking curve analysis and 
are summarized in the Table 2. X-ray standing wave analysis gave an additional information: the distance 
between the first (as counted from the substrate) atomic plane of the film and the diffracting plane of the 
substrate nearest to the interface (denoted as P hkl ) and the static Debye- Waller factor e~ w ° of the film. 

Comprehensive study of different type of garnet crystals with different kind of films on the surface were 
analysed with XSW method monitoring photoelectron yield and fluorescence radiation in |74|, [7^, [76| (first 
measurements of fluorescence radiation from garnet single crystals were performed in p3|]), Peculiarity of the 
garnet crystals is a comparably complicated unit cell containing 160 atoms (for the distribution of different 
atoms in the unit cell of garnet crystal in the (111) direction see Fig. 20(a)). Another interesting property 
is that due to the fact, that garnet crystals contain heavy atoms with Z » 1 propagation of x-rays in this 
crystals is highly reduced and fluorescence yield from this atoms is absorbed on a short distances of about 
1 -r 10/zm. 

The angular dependence of the total photoelectron yield from the gallium gadolinium garnet (GGG) 
crystals excited by the XSW field was measured for the first time in |7j (Fig. 20(b)). Due to a small 
value of the absorption depth Lq — smOs/fJ-o ~ 5.44/im, that is comparable with the extinction length 
L ex = in this crystal the reflection coefficient for (888) diffraction of CuK a radiation has a small 

~ 47% maximum value. At the same time the photoelectron yield curve has a phase sensitive dispersion like 
shape. 

In the same paper the angular dependence of the photoelectron yield for the GGG crystal with the 
epitaxial film (thickness of the film Ld = 2/jto) of the iron yttrium garnet (FYG) was measured (Fig. 
21). The same CuK a radiation and (888) reflection as in the previous case was used. For the theoretical 
description of the photoelectron yield in this case of thick film the general theory of bicrystal model described 
in the first part of this chapter can be applied. However for understanding the physics of formation of the 
photoelectron yield while the dynamical diffraction some simplified considerations can be used. Lattice 
parameters of the FYG film and of the GGG substrate does not differ essentially (the same is valid for their 
Bragg angles), however the Fourier components of the susceptibility Xh are quite different. For example, 
extinction depth for the FYG crystal is equal to L ex = 5.68 \xm that is essentially bigger then the thickness 
of the film (Ld << L ex ), so reflection from the film is kincmatical. In the kinematic approximation we have 
for the reflectivity, 



P R {0) 



y 2 (0) 



L, 



where angular deviation parameter y(9) is defined in ( |25| ) and C p 
factor. In this approximation we have for the photoelectron yield, 



(79) 

(1 + C 3 )/(l + C) is the polarization 



K p h{9) = 1 + K sin 2 



L„ 



y(0) 



y{0) 



y 2 {0) 



(80) 



Parameters in equations (|79|-pOD in the case of experiment [[74| has the following values K = 0.75, | e-r | = 
0.8 and we obtain for PR max = 0.092, K m i n = 0.75, K ma x = 1.33. This simple approximation fit quite well to 
the experiment. At the same time the photoelectron yield at the angular position of the substrate peak has 
the shape similar to the reflectivity curve (this is quite similar to the discussed before case of amorphous 
layer, however the physics of the process is different). Period of the standing wave formed in the substrate is 
different from that of the film (measured photoelectrons are excited only in the film). Due to this difference 
there is no correlation between the positions of the atomic planes on the escape depth of the electrons and 
the positions of the nodes and antinodes of the standing wave. In this experiment the escape depth of the 
electrons is about of L yi = 0.22 /im. On this depth due to the difference in the period of the film lattice 
and the substrate atom is shifted on 1.44A and at the same time the period of the standing wave for (888) 
reflection is equal to a\/3/24 = 0.89 A. Finally, due to a big variation of phase on the period of the standing 
wave field the third term in the expression for the photoelectron yield K p h(0) ( |77| ) cancel out and we obtain 
for K ph {6) = 1 + P R {6). 

In the same paper the fluorescence yield from the Gallium-Ncodymium Garnet (GNG) crystal and the 
films of FYG (thickness of the film L = 1.6/im) on the top of GGG crystal with XSW method were 
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investigated. Due to a big depth of yield of the fluorescence radiation {LP°- = 20. 4/xm for GaK a radiation 
and Ly^ = 15.1/xm for NdL a radiation) comparing to the extinction length (L ex = 1.54/im) the angular 
dependence of the fluorescence yield curves in the XSW field is determined by extinction effect discussed in 
the section 4.1. However as we can see from Fig. 22(a) the fluorescence yield from Nd and Ga atoms has 
different asymmetry on the tails of the curves. This is due to the fact, that for reflection (444) the phase of 
the standing wave on the neighbour atomic layers (on the distance of 1/24 of the period) differs by it. At the 
same time average coherent positions for Ga and Nd atoms coincide, but the effective coherent fractions F c 
are different: F^ a = 0.321 and F? d = 1.628. 

For the theoretical analysis of the fluorescence yield curves from the FYG films on the GGG crystal 
substrate the bicrystal model described in this section was used. While fitting in this case it was additionally 
taken into account the change of the susceptibility of the film Xfilm comparing to that of the substrate 
Xsubstrate as well as the change of the interplanar distance Ad/d. Results of the fitting are presented on Fig. 
22(b). In this case YK a fluorescence radiation from the atoms of the film and GaK a radiation from the 
atoms of the substrate in the case of (444) diffraction of AgK Q radiation were measured. Different shape of 
the curves reflects different conditions of the fluorescence yield formation in XSW field. Extinction depth is 
equal to L ex = 2.3pm and is bigger then the thickness of the film. The depth of yield of the YK a radiation 
is restricted by the thickness of the film and coincide with the volume of XSW formation in the film (the 
same as on Fig. 21, but with bigger value of Ld/L ex ). At the same time (when Bragg conditions for the 
film are fulfilled) the Ga fluorescence yield is decreased because x-rays being reflected in the film hardly 
penetrate to the substrate, where Ga atoms are located. Right maximum on the reflectivity curve on Fig. 
22(b) correspond to the reflection of x-rays from the substrate. In this case the angular dependence of the 
GaK a fluorescence yield has a big dip due to an extinction effect described before (due to a big escape depth 
of the fluorescence radiation comparing to an extinction depth). At the same time the yield of the YK a 
radiation in this angular region has a maximum of the same origin as described on Fig. 21. 

In the papers j7^, [7(| XSW method was used to study the positions of Bi 3+ ions in the lattice of a 
hcterocpitaxial film of yttrium-bismuth iron (Y 3 _ x Bi x Fe^Oi2) garnet. The main problem was to determine 
quantitatively the distribution of bismuth over the different dodccahcdral positions which are occupied by 
yttrium in a pure yttrium iron garnet. This information is essential for understanding the growth-induced 
magnetic anisotropy of this garnet crystals. The stated problem was solved by detailed analysis of the 
angular dependence of the Bi 3+ fluorescence. 

Ideal garnets have cubic symmetry with cations entering octahedral (a), tetrahcdral (d) and dodccahcdral 
(c) sublattices. In the Y-$- x Bi x F e^On system Bi 3+ and Y 3+ ions enter only the c sites. For the [001] growth 
direction there are two incquivalcnt groups of dodecahedral sites, the first consists of 16 sites, which are 
denoted as c xy , the second consists of 8 sites denoted as c z j77j. The distribution of c xy and c z sites within 
the elementary cell for this growth direction together with the distribution of the nodes and antinodes of 
(004) X-ray standing wave is shown in Fig. 23. Note that the sites belonging to the c xy and c z groups lie in 
different layers and consequently they can be, in principle, distinguished by the XSW method. 

For the analysis of the fluorescence yield from this garnet samples a bicrystal model consisting of the thin 
film and the substrate was used. If the depth of yield of the fluorescence radiation L y i is much smaller than 
the extinction length L ex the angular dependence of the intensity of fluorescence radiation in multicomponcnt 
crystal is given by Eq. (|7^) . For the reflection (004) the inequivalent groups of sites c z and c xy have different 
structure factors Sh(c z ) = —1 (8 sites), Sh{c xy ) = +1 (16 sites). The structure factor of the Bi 3+ ions, which 
enter the dodecahedral sublattice in Y->,- x Bi x Fe^O\i garnet films, is influenced by the distribution of these 
ions between c z and c xy positions. The fraction of Bi 3+ ions in c z positions will be denoted by parameter p 

Nf l 

P ~ N Bi T KfBi ' 

where N?\ N*j are the numbers of Bi + ions in c z and c xy sites, respectively. Then the structure factor 
of Bi 3+ ions for the (004) reflection is = 1 — 2p and due to Eq. (|7^) the fluorescence yield of Bi 3+ ions 
will depend on the value of the parameter p. The uniform distribution of Bi 3+ ions between c z and c xy sites 
corresponds to p = 1/3. The calculated angular curves of the reflectivity and _Bi 3+ -ion fluorescence yield for 
values of p in the interval to 1 are shown in Fig. 24. The angular yield curves of the Bi 3+ fluorescence 
demonstrate, as expected, a strong dependence on the distribution of these ions. For values of the parameter 
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p = (all Bi 3+ ions in c xy sites) and p = 1 (all £>i 3+ ions in c z sites) the amplitude of the structure factor 
of the Bi 3+ ions attains its maximum ( Sjf * = 1) and the phase ip 1 ^ 1 °f the structure factor equals and 
7T, respectively. The angular yield curves for these limiting cases show reversed positions for maxima and 
minima of the yield. When \S^ l \ decreases the ratio of the fluorescence to the background becomes smaller. 
Note that on the curve corresponding to p = 0.33 (uniform distribution of Bi 3+ ions over c z and c xy sites) 
a clear maximum occurs at the high-angle side (with respect to the maximum of Bragg diffraction on the 
film). At the same time the reflectivity curves (Fig. 24(a)) has very small dependence on the value of the 
parameter p. 

The diffraction curve for the sample studied is shown in Fig. 25. The angular distance between two 
peaks (corresponding to the substrate and the film) allows a direct determination of the lattice mismatch 
between the film and substrate Ad/d = 2.5 x 10~ 3 . In Fig. 26 the experimental angular dependencies of the 
fluorescence of Bi 3 + (Fig. 26(a)), Y 3+ (Fig. 26(b)), Fe 3+ (Fig. 26(c)) and Gd 3+ (Fig. 26(c)) are shown 
together with the reflectivity curve (Fig. 26(e)) in the angular range corresponding to the diffraction on the 
film. Results can be understood on the basis of the theory described in the first part of this Chapter. Three 
main factors influence the shape of the fluorescence yield curve, namely the amplitude and the phase of the 
structure factor of the sublattice of the ions under study and the extinction effect for a given thickness of 
the film. The experimentally observed shape of the curve is determined by the fact that the Bi 3+ sublattice 
structure factor is close to zero and the extinction effect for a given thickness of film leads to the decrease 
of intensity. The best fit (shown on Fig. 26 by solid curve) was obtained for the following values of the 
parameter p = 0.44 ± 0.02, thickness of the film = 3/j,m, and Debye- Waller factors of the film and 
the substrate exp(— Wf) = exp(— W s ) = 0.93. Small mosaicity of the sample is also taken into account by 
additional convolution of the calculated curves with the Gaussian function with a = 3.8 arcsec. The obtained 
value of the parameter p represents the main result of the papers |75|, |7(| . Due to interference nature of the 
XSW method it gives the possibility to obtain the value of this parameter with high accuracy (see Fig. 26a), 
where Bi 3+ fluorescence yield curve with the value of p = 0.33 corresponding to the uniform distribution of 
Bi 3+ ions) is also shown for comparison. 

The same bicrystal model was used in the analysis of a recent experiment [ [78[ where the lattice constant 
difference of the isotopically pure 76 Ge grown as a film on the top of a natural Ge single crystal was measured 
as a function of the temperature. The variation of XSW in the film was monitored with the photoclcctron 
detector. 
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6 XSW in Laue geometry 



Up to now we were considering formation of Standing Waves in crystal and secondary radiation yield for the 
Bragg geometry of diffraction. Registration of SR yield in Laue geometry has it's own peculiarities different 
from Bragg geometry of scattering and will be discussed in this Chapter. XSW in the Laue case with the 
registration of SR yield from the exit surface have been studied in several papers concerning lattice-atom 
fluorescence yield |7{J, external | ]80| a nd internal Jl8j photoeffect. A Laue-case interferometer was used for 
the location of chemisorbed atoms ]8l|] . In this Chapter we will present mainly theoretical and experimental 
results obtained in |8^, , where the formation of the fluorescence-yield angular curves in the Laue case and 
the possibilities of this geometry for the location of impurity atoms were studied. The advantage of the Laue 
geometry of diffraction (Fig. 19) is the possibility to study the position of impurity atoms distributed through 
the bulk of a crystal. In this case the sensitivity of XSW in the Bragg geometry decreases due to extinction 
effect (see discussion in Chapter 4 of this review). Moreover, in the Laue case it is easy to use different 
reflections (including asymmetrical ones) to study the impurity-atom positions in different crystallographic 
directions. 



6.1 Theory 

The general theory of SR yield described in Chapter 3 is valid for the Laue geometry as well. The yield of 
the SR is determined by the expression similar to ( |52| ) 

L 

* a {0) = J(L/ dz W {l^oMI 2 + \E h (z,6)\ 2 + 
o 

+2CF?(z)Re [E%(z, 6)E h {z, 9) cxp[^ Q (z)]]} , (81) 

where 

FZ{z) = e- wa ^\e±\ (82) 

is coherent fraction and L is the crystal thickness. The type of the SR and the type of atoms emitting it 
are characterized by index a. The phase <p a (z) = 2irP® + hu(z) = ipj- + hu(z) and the factor exp[— VF a (z)] 
correspond only to the atoms of type a, (p2- = 2nP^ is the phase of the complex value Sj- (^) and P c a is 
coherent position. The function Pyi{z) describes the yield of SR emitted by atoms of type a located at depth 
z. In the case of Laue geometry and the fluorescence yield from the exit surface Pyi(z) = c a exp[/^^ (z — L)\, 

where = l/OjfiLyl), Ly\ is the depth of the fluorescence yield, c a is the concentration of the impurity 
atoms (see Eq. (J55J)) . 

The field amplitudes Eq s (z,9) and Eh s (z,9) for the crystal with the deformation field u(z) and Laue 
geometry of diffraction (in this case parameter 7/j > 0) can be obtained from the general TT equations (fTs)). 
In this Chapter it will be discussed in details the situation, when crystal can be considered as a set of two 
layers (bicrystal model) with different lattice parameters, perfection and atomic compositions. A crystal 
lattice of each layer is characterized by a constant value of the parameters Ad/d, exp(— W) and Xo,hi which 
describe the change of the plane spacing, decrease of the coherent scattering amplitude and the composition, 
respectively. 

Let us consider like in the previous Chapter a layer of thickness at the depth zq < z < zq + Ld- The 
phase (f(z) in ( |l8| ) depends linearly on z because Ad/d is a constant within a layer 

ip{z) = ip{z ) + 2y v> Az/L ex , (83) 

where Az = z — zo, y v and L ex are constant, L ex is an extinction length averaged over the crystal bulk. 
The solution of TT equations (|l8| ) should satisfy the given values of Eo s (z, 9) and Eh s (z, 9) on the entrance 
boundary of the layer. 

As it was shown in Chapter 2, if the amplitude function R(z,9) is defined in the form (|23|), then it 
satisfies non-linear equation (|2^) (for Laue case lower signs has to be taken) with a boundary condition 
R = Rq on the upper boundary of a layer at z = zq. The solution of Eq. ( |24| ) in the Laue geometry with 
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constant parameters can be found analytically (see Eq. (114) from Appendix, where in the case of one 
layer the value of R n on the upper boundary has to be changed by Rq). Intensity of the transmitted beam 
Iq(z, 6) = \Eq(z,0)\ on the depth z defined in ( |34| ) also can be found analytically (see Eqs. (12C-123) from 
Appendix) . 

The contribution of the layer with the thickness Ld to the SR yield in the frame of these model is 



analysed in details in Appendix and is summarized in Eqs. (129 - 13l| ), where thickness of the layer d n has 
to be changed by Ld- This equations have to be complemented by the boundary conditions on the entrance 
surface of crystal, where Eq = 1 and R = 0. X-ray transmission (Pt(0)) and reflectivity (Pr(8)) are given 

by 



Pt{6) = \E (L,6)\%P R (6) = P t (6)\YR(L, 



(84) 



Finally the problem of the SR yield from the exit surface of the crystal in Laue geometry is completely 
defined. If measurements arc made in a double-crystal scheme with the first crystal-monochromator convo- 
lution with the first crystal reflectivity curve (according to (p^) has to be taken into account. 



6.2 Experiment 

Theory of the SR yield in Laue geometry of diffraction described in the previous subsection was applied in 
[ p2| . p3[ for the analysis of the experimental results. 

In the first series of measurements Si single crystals of different thickness with (100) surface orientation 



uniformly doped during growth with germanium (Nq c = 7.5 x 10 
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3 ) were used. GeK a fluorescence 



from the exit surface of the silicon crystals in the Laue geometry and (111) reflection was measured by an 
energy-dispersive Si(Li) detector (see Fig. 27). The angular dependence of GeK a fluorescence yield from a 
thick crystal with thickness 2.2 mm (Fig. 28a) shows a large maximum slightly shifted to the low-angle side 
with respect to the maximum of the diffracted intensity. On the contrary, in the case of a thin crystal with 
thickness 0.49 mm, the fluorescence curve (Fig. 28b) shows a dip near the center of the reflection range and 
a weak maximum at the high- angle side. 

According to our previous theoretical approach in this case the whole crystal can be considered as one 
layer with thickness Ld = L and fluorescence yield can be determined by Eq. ( |l29| ). In experiments the 



following conditions were also fulfilled Lq >> L y i >> L ex , where Lq = 7o//^o- Since L 



J IJ< 



>> L exi the third 



(interference) term in (129) is very small (IV'sl << \ipi\ \ip2\) and can be neglected. On the other hand, since 
Lq » Lyi the normalized fluorescence yield does not depend on L ex and L y i. Taking this considerations 
into account from (|129]) an approximate expression can be obtained 



k(L,0)~ ^ 4(0) exp [A/i fc (0)L] 



(85) 



fc=1.2 



where 
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2C^fpe a h exp(-W a ) cos(27rP c a ) 



(86) 



To derive (|S5|) the following values of the parameters were used Y = 1, y v = and expansion for the complex 
parameters x k in the power series of a small value Xi/Xr were used. 

The most significant feature of Laue-case X-ray diffraction is that two types of standing waves are 
formed in a crystal (Fig. 27b). One of these waves is a weakly absorbed field with the nodes on the atomic 
planes, corresponding to the well known Borrmann effect. The other one is a strongly absorbed field with 
the antinodes on the atomic planes. So, as follows from (|85|), the fluorescence yield angular curve differs 
distinguishably for a thick crystal (L >> Lq) and a thin one (L < Lq). The experimental curve shown 
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in Fig. 28a corresponds to L = 4.3Lq. In this case the SR yield is excited only by the weakly absorbed 
X-ray standing wave field [k = 1 in (|85|)]. The other feature of the Laue case is that with rocking a crystal 
through the reflection position, the standing wave fields do not move with respect to the atomic planes (as 
in the Bragg case), but only intensities of these fields increase or decrease. The weakly absorbed field has 
a maximum intensity at 9 = 9b- Even for impurity atoms lying strictly in the crystal nodes (P" = 0) one 
can observe an increase of the fluorescence yield at 9 ~ 9b in comparison with a background yield. At any 
displacement from the crystal node, the impurity atom occurs in the region of the increased field intensity 
and fluorescence yield increases even more. 

The obtained experimental curve (Fig. 20a) is in a good agreement with the theoretical one for the model 
of substitutional impurity atom position. The maximum normalized fluorescence yield is equal to 3.3. For the 
impurity atom between the reflecting planes (P" = 0.5) this value will be 13.9. Note, that the sensitivity to 
the position of the impurity atom increases with increasing crystal thickness. Calculated dependencies of the 
normalized fluorescence yield maximum on the impurity position with respect to (111) and (022) reflecting 
planes arc shown in Fig. 29 for different crystal thickness. Crystals with intermediate thicknesses are more 
appropriate for impurity location studies, due to a rapid decrease of the background fluorescence yield with 
increasing of crystal thickness. It should be noted that a randomly distributed fraction of impurities also leads 
to an increase of the fluorescence yield. Indeed, the static Dcbye- Waller factor for impurity atoms exp(— W a ) 
and the term cos(27rP") describing coherent displacements are included in ( ]85| ) as the multipliers. 

For a thin crystal, because of the excitation of both standing wave fields the situation is more complicated. 
Both weakly (k = 1) and strongly (k = 2) absorbed fields can make a significant contribution to the 
fluorescence yield curve, but at different angular positions: field 1 at 9 < 9b and field 2 at 8 > 9 B - The 
main factor now is a degree of interaction of impurity atoms with the standing wave fields but not the 
anomalous X-ray transmission. If the impurity atom is located on diffraction plane, it interacts with field 
2 more than with field 1. So, the maximum of the yield will be observed at 9 > 9b (curve 1 in Fig. 28b). 
The minimum on curve 1 occurs because field 1 does not interact with impurities and field 2 is not excited 
in this angular region. If impurity atoms are randomly distributed [in (|86j), exp(— W a ) = B = 0], this effect 
is compensated entirely by increasing the interaction with field 1 (curve 2 in Fig. 28b). In Fig. 30, one 
can see the fluorescence yield curves calculated for different positions (P c a ) of impurity atoms. It is obvious 
that with displacement of the impurity from the diffraction plane interaction with the weakly absorbed field 
increases and with the strongly absorbed field decreases. So, the fluorescence yield increases at 9 < 9b and 
decreases at 9 > 9b- The experimental curves both for thick and thin crystals unambiguously show that 
germanium is substitutional in silicon. Such behavior corresponds to the isovalent nature of this impurity. 

In the second series of measurements Si crystals with thickness 0.35 mm, (111) surface orientation with 
an epitaxial layer with thickness 1.6 /im on the exit surface were studied. During growth the epilayer was 
doped with boron and simultaneously with germanium {Noe = 10 28 rn~ 3 ). The experimental curve of GeK a 
fluorescence yield from the epilayer in the angular range of diffraction in the substrate is shown in Fig. 31a. 
In this case fluorescence was excited by XSW formed in the bulk of the crystal. An X-ray rocking curve in a 
wide angular range is shown in Fig. 31b. From this rocking curve the deformation in the epilayer comparing 
to that in the substrate can be estimated as Ad/d = —3.8 x 10 -3 . In the conditions of X-ray diffraction on 
the epilayer, the so-called 'kinematical' standing waves are formed in this layer. The fluorescence yield in this 
angular range is shown in Fig. 31b. Weak minimum and maximum were clearly observed on the increasing 
background. The angular dependence of the background is due to the difference between the plane wave and 
the wave incident on the layer due to the diffraction in the bulk. 

In this series of experiments the layer with thickness Ld = 1.6 fj,m is much less than the extinction length 
L ex = 12.6 fim and L > Lq. Impurity atoms which are only in the layer 'see' both standing wave fields 
and their interference in the angular range of the X-ray diffraction in the substrate. The interference term 
oscillates with increasing z and the period of these oscillations is different at different values of the incident 
angle: at 9 = 9b the period is equal to nL ex and it decreases with increasing \9 — 9b\- Since L » L ex , the 
fluorescence yield angular curve for the case of the incident plane wave strongly oscillates with the period less 
than the angular divergence of the real incident beam. So, in the treatment of the experimental curve the 
interference term can be averaged again giving a zero contribution, although the reason for this averaging is 
quite different for that of the previous case. 

Taking all this into account, it is clear why the calculated curves in Fig. 31a for substitutional and 
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randomly distributed impurities are very similar to the curves in Fig. 28b, although the calculations were 
carried out in a double-layer model including convolution. In this case opposite to the previous one the 
experimental curve is in a good agreement with the model of random distribution. This means that impurity 
atoms are in all possible positions with respect to the standing waves formed in the substrate. Nevertheless, 
impurities in the layer could be ordered, but with period different from the period of the substrate reflection 
planes. 

The problem of the impurity position can be solved by measuring the fluorescence yield curve in the 
angular range of X-ray diffraction on the epilayer. Since the reflection of X-rays by the epilayer is small 
(L d << L ex ), the standing wave formed in the layer has the intensity weakly oscillating around the average 
value. It is the case of so-called kinematical standing waves. The total reflectivity amplitude R(z, 9) (^3|) in 
the kinematical approximation can be presented in the following form 

R(z,9) = Ce -«' e .^ + ( 87 ) 
R (A9 + A9 B )e 2iqz , 

where the coordinate z starts from the upper boundary, q = y/L ex = (n sin 29b / \"fh)A9, A9 = 9 — 9b is 
the angular deviation from the Bragg angle for the layer, A9b = 9b — 9 b , where 9b and 9 B are the Bragg 
angles for the layer and the substrate, i?o is the reflectivity amplitude of the substrate. 

The fluorescence yield from the layer is defined by the general equation (^Tj) but at L y i >> L d one can 
put Pyi{z) = 1 in (fjl]). Moreover, one can consider that the amplitude of the 'incident' wave docs not 
depend on z and this is the component of one of the standing waves formed in the substrate. Since X-ray 
diffraction by the layer is observed at 9 > 9 B the 'incident' wave is the component of the strongly absorbed 
field (k = 2). In the case of the paper [^3) the shift of the Bragg angle in the layer with respect to the 
substrate A9b is large enough. It allows to neglect the h— component of the substrate standing wave and 
put Ro(9) = in the following consideration. 

For simplicity two different cases will be discussed: the first one is for randomly distributed impurities 
(F c a = 0) and the second one is for impurities in crystal nodes (P" = 0). 



In the first case, substituting (|87|) into ( 81 ) and performing integration the following result is obtained 



K (L d ,9)~ \E (9)\ 2 L d 



1 | P fCe- w \ 2 / 



2 \ qLex J V 2 1 L d 



As follows from (p9|), on the background yield which increases with increase of angle 9 in the range 9 > 9 B , 
a very weak maximum should be observed with symmetrical shape and intensity proportional to (L d / L ex ) 2 . 
In the second case, 



K{L d ,A9) ~ \E (9)\ 2 L d 



-w 

(89) 



1 + 2^F-Ce~ w / _ sin(2 g £ d ) 

qL ex V 2qL d 



From (|84|) quite different angular dependence of the fluorescence yield follows. It is determined mainly by 
the interference term and after substraction the background it has a dispersion form with the maximum at 
A9 > and the minimum at A9 < both being proportional to (L d /L ex ). Just this behavior is observed 
experimentally (see Fig. 31b) pointing to the correct positions of germanium in the layer lattice. 

The theoretical curves shown in Fig. 31 were calculated using the general approach described in the 
beginning of this Section. Both experimental fluorescence yield curve and X-ray rocking curve are weaker 
and broader. This effect can be explained by formation of the mosaic structure due to a high concentration 
of germanium in silicon. 

Finally, in this chapter we have demonstrated, that standing waves giving useful structural information 
can be formed in thin and disordered crystals. In the most of the previous papers it was considered that 
XSW fields can be formed only in large and nearly perfect crystals. 
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7 Model of a multilayer crystal 



7.1 Theory 



The previous case of a bicrystal model in Bragg and Laue geometry of scattering is the simplest (however 
quite useful for the analysis in many situations as it was demonstrated in the previous Sections) model of a 
crystal with deformed surface layer. At the same time there are many cases, when the profile of deformation 
and amorphization has a more general functional form then just a step function. In this case different 
approaches can be used to solve the problem of SR yield from such a crystal. If the profile of deformation 
and amorphization is known, then one can solve TT equations either in the form of a set of equations ( |l8| ) 
or in the form of Eq. ( p4| ) for the amplitude R(z, 9) numerically and then substitute the obtained solution 
into Eq. ( p2| ) for the angular dependence of the SR yield k{9). Such approach was proposed for the first time 
in the paper for the analysis of the photoelectron yield from the crystal with the deformed surface layer 
(see Fig. 8). Later the same approach was used in |7l| for the analysis of the wave fields in a crystal with 
the deformed surface layer taken in the form of a step Fermi function. 

However, in many cases another approach [ p4[ can be useful. One can divide deformed (and amorphisized) 
layer with an arbitrary shape on a several sublayers (see Fig. 32) with a constant structure parameters in 
each layer (these are in our case deformation Ad/d and amorphization e~ w ) and to solve TT equations in 
each sublayer analytically. Then substituting this solutions in the general equation (^2|) for the SR yield the 
desired result for a multilayer crystal can be obtained. 

The details of this approach are given in the Appendix. Analytical solution of Eq. (^J) for the amplitude 
R(z, 9) in the layer z n < z < z n + d n in the frame of the described model can be presented in the form of 



Eq. (114). In this equation both geometries of diffraction (Bragg and Laue) are considered simultaneously 
From the obtained solution it is easy to get standard recurrent equations. Defining by Ri n the reflection 
amplitude for both cases on the "known" boundary (Ri n = Rd n for Bragg case and Ri n = R n for Laue), 
by Rout the reflection amplitude on the unknown boundary and by r out solution for Ri n = (reflection 
amplitude by thin layer) we have from (114- 117) 
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where T is the transmission amplitude. 

Directly from Eq. ([52|) solution for the direct beam Eq(z, 9) in the layer z n 
in the form 

E o (z,0) = T(z,e)E o (z o ,0), 

where T(z,6) is the transmission amplitude through the layer defined in 
the layer 1(9) = \E a {d, 9)\ 2 is determined by Eq. ( p^) . 

According to (|5^) the normalized SR yield from the n-th layer in general case of multiatomic unit cell 
for the sort of atoms a and with the exponential form of the yield function (|5^) from the entrance surface 
can be presented in the following form 



121 



(92) 

Intensity transfer through 



K a (9) = xto J dzP yl (z)\T(z,9)\ 2 x 



l + m?4& + 2CRe 



\E (z,9)\ 2 

Substituting to this expression previously defined functions we obtain 

d 

k%9) = x-o J ' dzP yl (z)\T(z,9)\ 2 x 
o 



Eh(z, 0) i V °-( z )-w a 
E (z,9) h 



(93) 
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In this equation two important parameters of the theory i.e. coherent fraction (F£) and coherent position 
(P c a ) for the atoms of the sort a in the n-th layer are introduced. They are defined as following 



= \^W- W , <p a c = = A^ a + arg( £ |) (95) 
and have a constant value in each layer, but can change from layer to layer. 



Taking into account analytical solutions for the amplitudes R(z,9) (114). T(z,9) (121) and exponential 
form for the probability yield function P y i{z) ( p5| ) integral in expression for n a {9) can be taken analytically 



(see Eq. (129) in Appendix). Finally the total yield of the SR from the whole multilayer crystal is the sum 
of the yields from every layer 

N 

< o4 (0) = $>„(0K(0). (96) 

71=1 

where z\ — 1 and z n (9) = z n ^i(9)\T(d, 9)\ 2 e~ flyid . So, the whole algorithm for calculation of the angular 
dependence of the SR yield from the multilayer crystal in both (Bragg and Lauc) geometries is formulated. 
In the following several examples of application of this approach will be given. We have to note here, that 
described in the previous Sections bicrystal model is a particular case of this more general approach with a 
number of layers N = 2, where the second (thick) layer can be considered as a substrate and the top (thin 
layer) as an epilayer. 

7.2 Applications 

Described above approach of multilayer crystal was used in the papers |85| , |S6| where silicon crystals were 
implanted with Fe and Ni ions and analysed with XSW technique and Rutherford backscattcring (RBS) 
method. The XSW technique was applied previously [|[ ^8, 87 for structure investigation of Si crystals after 



implantation with the different impurity atoms (As, Bi). However, in this papers more simplified model of 
the damaged crystal was assumed. 

Below, following the papers |8^, ^6|, we will discuss in details results of XSW analysis of silicon single 
crystals implanted with Fe ions (dose 5 x 10 15 ions/cm 2 at an energy of 80 keV). After ion implantation 
samples were annealed at 750° in the N% atmosphere. The fluorescence yield from Fe atoms was measured. 
Results of this measurements on the samples before and after annealing are presented on Fig. 33. 

The RBS measurements were performed with 2 MeV 4 He ions. On Fig. 34 the results of the damage 
analysis from the RBS spectra are presented. On this figure the evolution of the Fe atoms depth distribution 
is also shown. The measure of the damage density k is given by the ratio of the displaced atom density to 
the crystal atom density , hence k — 100% denotes a complete amorphization. An important feature of the 
Fe— atom depth distribution for an as implanted sample is its quasi- Gaussian shape with the maximum (i.e. 
Fe mean projected range) located at the depth of 70 nm and range struggling of 80 nm. As it follows from 
the analysis the region of amorphization extends over 190 nm. After annealing the Fe distribution become 
narrower, however, its maximum position did not change. Also the peak area has remained the same. This 
effect is apparently due to the low solubility of iron in silicon. In such a case the regrowth front pushes away 
impurity atoms in concentrations exceeding the solubility limit. Due to this process the impurity segregation 
occurs in a not yet regrown region. 

For the theoretical fitting of the angular dependencies of the fluorescence yield for implanted region of 
the silicon sample a "layered" model described in the first part of this Chapter was used. In the frame of 
this model each deformed layer was characterized by its constant value of thickness d, deformation Ad/d 
and amorphization factor exp(— W). Distribution of foreign atoms in these layers was characterized by two 
main parameters: Pf e — the coherent position, i.e. the position of the mean plane of the impurity atoms 
with respect to the diffraction planes of the crystal and Ff e - coherent fraction, which describes the static 
and thermal displacements of the atoms from the mean position P^ e (see for definition Eq.(p5|)). Naturally, 
the contribution of each layer to the fluorescence yield is determined also by the concentration distribution 
n^ e of the impurity atoms in the surface layer of the silicon sample. All other parameters were the same in 
all layers. 

For the theoretical fitting of the fluorescence yield from the implanted crystal before annealing 14-th 
layer model and after annealing 11-th layer model were used. The profile of amorphization e~ w of the 
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silicon surface and the distribution of implanted Fe atoms nf e for each layer were taken from the results of 
the RBS experiment (Fig. 34). The best fitting in the case of not annealed sample was obtained with the 
value of the coherent fraction Ff e = 0.00 ± 0.02. This result, as it was expected, mean that the implanted 
atoms occupy random positions in the lattice cell. Parameters of the layers and results of the fitting for 
the sample after annealing are summarized in Table 3. In addition to results obtained from RBS it was 
proposed that some parts of implanted atoms after annealing locate in substitutional positions of Si atoms 
(parameter Pf e = 0). During the fitting procedure the coherent fraction of Fe in the corresponding layers 
Ff e and the value of deformation Ad/d of the surface layer were varied. To minimize the number of fitting 
parameters the uniform distribution of Ff e and Ad/d in the surface layer for the depths ~ 40 -j- 120 nm was 
taken. The best fitting result (see Fig. 33 and Table 3) gives for the average value of the surface relaxation 
Ad/d = (-9.6 ± 0.09) x 10~ 4 and the value of the coherent fraction F c Fe = 0.27 ± 0.02. This value of the 
coherent fraction means that approximately 27% (with respect to the diffraction planes) of implanted Fe 
atoms are in substitutional Si positions. It is interesting to point here that from the independent XSW 
analysis it was obtained the same value of substitutional fraction of Fe atoms as after the RBS analysis. 

The obtained results show the effectiveness and mutual complementarity of the two methods for investi- 
gation of implanted crystals. The RBS experiments directly give the depth profiles of damage density for Si 
surface layer and the profile of implanted ions concentration. The XSW analysis is the precise phase sensitive 
method. It gives directly the value of the phase <p(z) that in our case is equal to <p(z) = hu(z) + <p£ e , where 
u(z) is the deformation of the surface layer (in the case of the uniform deformation we have u z = (Ad/d)L 
and iff e = 2irPf e . As one can see from this expression the phase is the sum of deformation and the position 
of the impurity atoms in the unit cell and additional information is necessary to determine this two values 
independently. As it was shown in our analysis the combination of XSW and RBS methods gives us the 
opportunity to solve, in principle, this problem uniquely and determine positions of impurity atoms in the 
unit cell as well as the total deformation (relaxation) of the surface layer. 

In the recent paper (88) for the analysis of the fluorescence yield from the superlattice crystal was also 
used a similar multilayer model. 
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8 Crystals with an extended deformation field 



Up to now we were mainly considering the cases when the deformation field is located in the thin subsurface 
layer. However this condition is not always fulfilled. There exist a big class of deformations in crystals, that 
extends through the whole volume of the crystal. These are, for example, crystals with the Uniform Strain 
Gradient (USG), bent crystals, vibrating crystals and etc. In the following of this section we will show how 
x-ray standing wave analysis can be performed in such crystals. 



8.1 Crystals with the uniform strain gradient. Bent crystals 

Theory of x-ray standing wave behavior in crystals with USG and as a special case bent crystals was 
for the first time developed in papers [^9], 90, Later similar theoretical approach for the analysis of 
the crystals with USG was performed in [92|. Investigation and characterization of crystals with USG 
(bent crystals) by different diffraction methods is of great practical interest - first of all, due to their wide 
application as x-ray monochromators at synchrotron radiation sources [B3l H| or focusing spectrometers for 



x-ray plasma diagnostics 95 1, and secondly, due to the fact that a lot of typical samples (e.g. implanted 



crystals, heterostructures) are often curved. Below we will present results obtained in |8^, |9(J, 91. 

The problem of the x-ray standing wave analysis in bent crystals is quite complicated due to the fact 
that the deformation field in such crystals is in general a two-dimensional function and extend through the 
whole volume of the crystal. Below we will consider a special case of such field: crystal with USG with the 
deformation field hu = Bz , where z is the depth from the crystal surface. If the incident wave is spherical 
and the source is on the Rowland sphere then cylindrically bent crystals arc particular case of crystals with 
USG 

The general problem of two variables in this specific case is reduced to determination of the amplitudes 
Eo(z) and E/j(z), which depend on one variable z. The x-ray amplitude R{z,9) defined, as before, in the 
form (|2^) satisfies equation (|2^). In the case of the crystal with the USG an effective angular deviation 
parameter y e fj(9, z) is a linear function from the depth z: 

y eff (6, z) = + y v {z) = y{6) + C{z/L ex ), (97) 



where y(0), y<p(z) are determined by the incident angle on the surface of the crystal (25) and the deformation 
field (|2^) and C is the dimensionless parameter defining the strength of this field. For the cylindrically bent 
crystal the relationship between parameter C and the curvature radius R is given by p9] 

where K is an elastic constant. 

For the analysis of the angular dependence of the SR yield in the crystal with USG general equation ( |52] ) 
can be used. The type of the secondary process (fluorescence, photocmission) will be determined, as before, 
by the imaginary part of susceptibility \oi an< i probability yield function P^(z). 

Analytical solutions for the wave fields Eo(z) and E^(z) in the crystal with USG were obtained in 
|p6| , however they can be effectively analysed only in their asymptotic form. The problem of x-ray Bragg 
diffraction and standing wave determination in a bent crystal can be effectively solved with numerical 
methods. In the papers |89| , |90|] the modified layer model described in the previous Section was used (see 
Fig. 35). The whole deformed crystal was divided into layers, each of thickness Az = L ex /(mC), where m 
is an integer. In each layer the value of the angular deviation parameter j/j was considered constant and its 
value differs from layer to layer by the magnitude Ay = 1/m (parameter m was taken from the condition of 
small change of Ay from layer to layer). Analytical solution of Eq. (E4h for i?(z,y i ) in each such layer and, 



hence, in the whole crystal are given by (114). The angular dependence of the SR yield in the crystal with 
USG can be obtained straightforwardly by substituting the solutions for R{z,yi) into Eq.(^), integrating 
it in each thin layer with the exponential probability yield function and summing over all layers (see Eqs. 

(US))- 

The angular dependencies of the reflectivity and SR yield were calculated according to the described 
model for the different values of the parameter C. It was taken equal to 0, 0.1, 1 and 10, which, according to 
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(Pq), correspond to the curvature radius R equal to oo, 550 cm, 55 cm and 5.5 cm respectively. On Fig. 36 
results of these calculations in the specific case of (400) symmetrical Bragg diffraction, CuK a , a— polarized 
radiation in a thick Si crystal are presented. 

The behavior of the x-ray reflectivity (the broadening of the curves in the central part, oscillations on the 
tails of the curves) correspond completely to the behavior described for the first time in ^9|. The angular 
dependencies of SR yield were calculated for the processes with the different depths of yield. Two cases 
were considered: fluorescence yield of SiK a radiation (with L y i = 70 y,m) and photoelectron yield (with 
L y i = 0.1 /jra) from the atoms of the crystal matrix. 

The behavior of SR at small values of parameter C < 0.1 (Fig. 36b) in the range of the angular 
parameter \y\ < 1 differs only slightly from those for the perfect crystal (Fig. 36a). The photoemission 
yield {L y i « L ex , L ex = 3.6 [im) for this angles has a dispersion like behavior and is quite sensitive to the 
phase shift. On the contrary, for the fluorescence yield (L y i << L ex ) due to extinction effect in the same 
angular region the characteristic dip is observed, which suppresses the phase sensitivity. However, already at 
small values of curvature characteristic oscillations can be seen on the tails of SR curves. They correspond 
to oscillations on x-ray reflectivity curve. With increasing curvature one can sec first of all a broadening 
of the phase-sensitive region and secondly the increasing of the amplitude of the oscillations on the tails of 
the photoemission yield curves (Fig. 36c). For large values of curvature the SR yield curves become, on the 
average, flatter and flatter and the period of oscillations is increased (Fig. 36d). These curves correspond to 
a strongly deformed crystal. So, in the case of the curved crystal we have an interesting possibility to study 
on the same sample the standing wave behavior in both limits of a weakly and strongly deformed crystal. 

The x-ray standing wave field behavior in the bulk of the deformed crystal is determined first of all by the 
type of deformation in the crystal (see for e.g. |io|, |7^]). Reflectivity Pn{z,y) and x-ray standing wavefield 
in the bulk of the crystal with USG were analysed in details The intensity of the wavefield on the depth 
z has the form 

I(z,y) = \E (z,y)\ 2 \l + ^YR(z,y)e lA ^\ 2 (99) 

and the phase shift A^s" is determined by the positions of the impurity atoms relative to diffraction planes. 
The values of Eq(z, 8) and R(z, 6) were calculated according to the model described before. 

On Fig. 37a two-dimensional distribution of the reflectivity Pn(z,y) in the {z, y) plane is presented. 
Calculations were made for the parameter C = 1 and for the depths z = -f- 2L ex . As it follows from the 
obtained result for a crystal with USG the reflectivity at an arbitrary depth z (while remaining invariable 
in its shape) is shifted in the angular range according to (|97]). Such behavior of the reflectivity indicates 
that at the depth z there exist an angular region of the total reflection of x-rays where the phase of the field 
amplitude changes from to n (so-called phase-sensitive region). For the fixed incident parameter yo{9) that 
region is located at the depth z e ff ~ \yo{8)L ex /C\ and has an effective thickness Az e ff ~ \2L ex /C\. Thus 
the size of this region is determined by the value of the parameter C. 

The two-dimensional distribution of the intensity I(z, y) ( |99| ) in the {z, y) plane for the same value of the 
parameter C arc shown on Fig. 38. On Fig. 38a intensity of the wavefield is calculated on the diffraction 
planes (Ay>° = 0) and on Fig. 38b between the diffraction planes (Ay>° = w). Comparing these intensity 
distributions one can notice that they differ both by the value and by the position of oscillations the tails 
of the curves. This effect of phase sensitivity of the standing wave intensity (|99| ) in the bulk of deformed 
crystal (contrary to the case of perfect crystal) can be understood if the behavior of the transmitted wave 
intensity T 2 (z,y) = \Eo(z,y)\ 2 in the crystal with USG is analysed. 

The two-dimensional distribution of the transmitted wave intensity T 2 (z,y) in the (z,y) plane for the 
parameter C = 1 is presented on Fig. 37b. According to the analysis performed in Section 4 the phase 
sensitivity of the standing wave intensity in the bulk of the perfect crystal at depths z ~ L ex is suppressed 
by the extinction effect. The situation is quite different in the case of the crystals with USG (see Fig. 
38b). This more complicated picture is a result of interference of the incident and diffracted waves with 
a deformation field of a crystal. This interference gives rise to an increase of the transmitted wave on 
the edge of the phase-sensitive region even at the depths z ~ L ex (Fig. 37b). Just due to this effect the 
phase sensitivity of the standing wave intensity in the bulk of the crystal with USG is increased (Fig. 38). 
Therefore, although for the perfect crystals applications of the standing wave method arc limited to thin 
subsurface layers in the case of the crystals with USG this method can be extended to investigation of thin 
buried layers (monolayers) of impurity atoms even located at the depths z > L ex . 
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Different caicuiation approach was used in the paper |n] (see also (]92||)) for the analysis of x-ray standing 
wave field on the surface of the crystals with USG (deformation field was taken in the same form hu(r) = 
B ■ z 2 , where B is a constant). According to this approach the amplitudes of the x-ray field Eo(r) and 
E/j(r) were obtained from the direct integration of the TT equations in the form (|l^). For the details of this 
calculation algorithm in the Bragg geometry of diffraction see, for example, review paper ]97| . 

In the paper |9l| the angular dependencies of x-ray reflectivity and intensity of x-ray standing waves ( |99| ) 
on the surface of the crystal with USG were calculated for the case of (400) symmetrical Bragg diffraction, 
CuK ai o — polarized radiation in a thick crystal of silicon. Results of these calculations are presented on 
Fig. 39. The XSW curves on the surface of the crystal were calculated as before for two different values of 
the phase Aip^, = and Atp^ = it that correspond to coordinate r on the surface of the bent crystal, that 
either coincide with the position of the diffraction planes or to the position between diffraction planes. In 
the calculations the USG constant B was taken equal to B = 2 cos 2 8b O-/Lex)0c, where f3 c = tt/(2L^ x ) is the 
critical value used in the theory of x-ray propagation in bent crystals |97fl and L^ x is the extinction length 
in the symmetric Laue case. This value of constant B correspond to parameter C = 1 from Eq. (^) . The 
XSW curves in the deformed crystal comparing to that in the perfect one (Fig. 39) shows broadening in the 
central part (in the so-called phase sensitive region) and also a sharp oscillations on the tails that correspond 
to oscillations on the reflectivity curve. The angular position of these oscillations depend first of all on the 
strength of the deformation field (determined now by the value of parameter B) and also on the phase shift 
value Aip®. The origin of this oscillation is the waveguide character of X-ray propagation in crystals with 
USG p9j . It is important to note here, that the standing wave picture obtained from the direct numerical 
calculations of the TT equations for the crystals with the USG are in a good agreement with the simplified 
multilayer model used in [^9[ |9(| and described in the beginning of this Chapter. However, from the point 
of the waste of the computer time the first method is much more effective. 



8.2 Vibrating crystals 

Another case of x-ray standing waves in the crystals with extended deformation field: namely vibrating 
crystals was considered in the paper |9q| . If in the previous subsection the deformation field in crystal 
can be controlled by the strain gradient or radius of curvature in the case of vibrating crystal this can be 
achieved by changing the amplitude and the wave vector of ultrasonic excitations in the crystal. What is 
more important for x-ray standing wave technique - extinction length in vibrating crystals can be also tuned 
in controllable way. The periodic deviation of the reflection planes from their ideal position in a perfect 
crystal misorientates the inner crystalline structure of the crystal, effectively " switching off" the whole parts 
of the crystal from the dynamical diffraction scattering. Increasing of the ultrasonic amplitude increases the 
deformation field and due to its periodic structure creates new dynamical reflexes (satellites) with extremely 
big effective extinction length. In a limit of a big ultrasonic amplitude, one can achieve the kinematical limit 
of diffraction when the initial reflecting power of the crystal is distributed over the large angular interval. 

Below, following |98| , the angular dependence of the SR yield will be considered for the special case 
of vibrating crystals, crystals disturbed by longitudinal ultrasonic vibrations propagating normally to the 
crystal surface. The advantage of this scheme is the uniform distribution of the deformation field along the 
crystal surface ant its layered structure in the bulk of the crystal. 

In this case a vibrating crystal can be presented by a set of identical layers with thickness equal to 
the period of the ultrasonic wave. Every layer can be approximated by a combination of thinner layers 
with constant value of interplaner distance (see Fig. 40). This periodic structure give rise to an additional 
satellites in the reciprocal space with the position h + ?iK s , where K s is the wave vector of the ultrasonic 
wave. With increasing the amplitude of the ultrasonic wave the undisturbed part of the crystal (and the 
reflecting power of the main reflex) decreases whereas the fraction corresponding to the satellites increases. 
Below the special case of the so-called " high- frequency" ultrasonic wave, when the ultrasonic wavelength is 
smaller then the extinction length of the given reflection will be considered. This gives the possibility of 
investigating each satellite independently, because they do not overlap in the angular range with each other 
and with the main reflex. 

The displacement field of the ultrasonic standing wave will be taken in the following form 

u(z,t) = w cos cut cos K s z, (100) 
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where w and u> are the amplitude and frequency of the ultrasonic vibrations (the wave vector K s is directed 
along the normal to the crystal surface) . 

The dynamical symmetrical x-ray diffraction on a thick crystal with the diffraction vector h directed 
along the normal to the crystal surface (axis z) will be considered below. The x-ray wavefield amplitudes 
Eq(z 7 9) and Eh{z,0) in the displacement field ( |10C| ) can be effectively found from the set of TT equations 
( |l8| ) . In this case Fourier component of the susceptibility of the deformed crystal Xk can ^ e presented as an 
infinite sum 

oo 

Xh = Xhe lhu{z} = Xh J n Q™)i n e X p(™K s z), (101) 

n— — oo 



where J„(hw) is the Bcsscl function of real argument. Each term in the infinite sum (101) correspond to the 



satellite with the diffraction vector h + nK s . At small values of the ultrasonic amplitude w only two first 



order satellites (n = ±1) are essential in the expansion (101). The amplitudes Eq(z,9) and Eh(z,9) can be 
also expanded in the same way 

oo 

E g (z,6)= ]T E g , n (z,6)exp(inK s z),g = 0,h. (102) 



Substituting expression (102) into TT equations (|lj) it is possible to obtain a set of similar equations for 



the amplitudes Eq iTI (z, 9) and Eh. n (z, 9) with the following substitution Xh — * Xh.n and a-*a„, where 

a n = a + nK s —^- and Xh.n = i"J n (hw)x?i (103) 

TT 

are the new values of angular deviation parameter and susceptibility of the n-th satellite. After straight- 
forward calculations analytical solutions for the amplitudes of transmitted (only component Eq i q(z,9) is 
important) and diffracted Eh n (z, 9) wave in the vicinity of the n-th satellite can be obtained (see for details 
@). 

The x-ray wavefield in the vicinity of the n-th satellite has the form of a standing wave with the period 
27r/|h + nK s | 

E ~ £ ,oe' ikor + £h. n e' i(kh+,lKs)r , (104) 

It is important to note here that the period of the standing wave does not depend effectively on the order 
of the satellite, because the diffraction vector is much bigger than the wave vector of ultrasonic wave. The 
only parameter varying from satellite to satellite is the effective extinction length. In the angular region of 
the n-th satellite the effective extinction length is equal to 

L - - m- ,105 » 



where L ex is the extinction length in the undisturbed crystal. According to (105) an effective extinction 
length is a function of the ultrasonic amplitude w and can be tuned by the change of the voltage applied to 
ultrasonic transducer. 

The secondary radiation yield from the vibrating crystal can be obtained in the frame of the general 
theory discussed in Chapter 3. Substituting the obtained values of the wave fields E$${z,9) and Eh, n {z,9) 
into Eqs. (|§ |^) the angular dependence of the total secondary radiation yield N„ (9) in vibrating crystal 
on the n-th satellite can be presented in the form of the infinite sum 

oo 

N n (9)= ]T N, hm (9). (106) 

m— — oo 

Each term in this expansion can be calculated analytically if exponential form d53) for the probability 



106 is 



yield function P y i{z) is assumed. Finally it was obtained [B8l that every term N n ^ m (9) in expansion 
proportional to the relative contribution of the scattering potential J m (hw) and through this Bcsscl function 
depend on the value of the amplitude of the ultrasonic wave. 
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This theoretical approach was applied for the calculation of the photoelectron yield (with the depth of 
yield L y i = 0.2^m) from Si (111) crystal CuK a radiation. The ultrasonic wavelength A s was taken equal 
to 0.7/xm. In this case general equations can be simplified and in the limit of small ultrasonic amplitude the 
photoelectron yield near two first satellites (n = ±1) is equal to (higher order terms in expansion ( |106| ) can 
be neglected) 

Ni(6) ~ iVi,i(0) + N lt0 {9) + (107) 

At the same time contributions to the main reflex from two first-order satellites are completely compensated 
and the total SR yield can be approximated by the expression 

N (6) N 0i0 (0). (108) 

Results of numerical calculations of the angular dependence of the photoelectron yield for the different 
values of ultrasonic amplitude hw = 0.1, 0.5, 0.9 are presented on Fig. 41. As it follows from the numerical 
calculations the photoelectron yield near the main reflex varies very slowly, however, the yield excited in the 
angular range of the satellites is very sensitive to the value of the ultrasound amplitude. 

It is important to note here that approach discussed in this subsection gives the intensity of the diffracted 
x-ray waves and the intensity of the SR yield for a fix time. In ordinary x-ray scattering experiment on 
vibrating crystals only the quantities averaged over the period of vibrations are measured. However results 
discussed here can be, in principle, observed in stroboscopic experiments, that are available on synchrotron 
radiation sources. The coincidence between the time-periodic flashes of the synchrotron radiation beam and 
the period of the ultrasonic vibrations gives the possibility to scatter the x-ray wave on the "frozen" distorted 
lattice. 

In the recent paper |)9] developed above theory was extended for the case of the yield of the thermal 
diffuse scattering in vibrating crystals. 
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9 Phase problem 



It is well known, that in any experiment based on the scattering of radiation by condensed matter, the 
measured quantity is an intensity, from which one can determine only the modulus of the complex valued 
scattering amplitude. In order to obtain comprehensive information about the structure of the scattering 
object, however, it is also necessary to know the phase of this amplitude. This is a general problem of 
scattering that arises in many applications, for example, in the study of the structure of the surface layers 
of single crystals by means of x-ray scattering. In this case the structure of the surface layer is directly 
connected via Fourier transform with the angular dependence of the scattering amplitude. At the same 
time a standard x-ray diffraction experiment does not give an opportunity to determine the phase of this 
amplitude. 

To solve this problem of phase determination indirect methods are widely used in x-ray crystallography 
(see for review of different methods a book [[lOOf ). During last years it was a rapid development of the 
direct methods for the phase determination of the structure factors, for example, by means of multibcam 
diffraction |38| [3E], 101 [. Another direct method of the structure analysis and phase determination is 



based on XSW method. As it was demonstrated in Chapter 4 coherent fraction ^(h) and coherent position 
P" (h) represent, in fact, the h— th amplitude and phase of the Fourier decomposition of the distribution of 



atoms under consideration (see also 58 611] ). Due to availability of the synchrotron radiation sources of the 



third generation new methods based on the utilization of the high coherency properties of the synchrotron 



radiation for the phase determination has become possible [102, 103 . In this case powerful iterative methods 



for the phase retrieval are used 104] 



Here we will discuss another approach for the direct phase determination of the scattering amplitude 
based on the simultaneous measurements of the reflectivity and SR yield in the large angular range. This 
information can be used for the direct and what is very important unique determination of the profile of the 
structural changes of the crystal surface layer. The idea of this method was for the first time proposed in 
p2| and later realized practically in |105| . 

If the thickness of the deformed layer in the sample satisfies the condition < < L ex , where L ex is 
the extinction length, then the scattering of x-rays by this layer is kinematic and the scattering amplitude 
R{9) is given by 

oo 

R{9) = \R{6)\ e la ^ =RoJ dzF{z)e lq{0)z , F{z) = /(z)e^ (z) . (109) 

o 

Here Rq = — z7rxh/(Asin0g), q(9) = (47r/A) cos^A^, f(z) — cxp(— W(z)) is the static Debye- Waller factor 



and ip(z) = hu(z). According to expression (10E), outside the region of the total reflection the scattering 
amplitude R(9) is the Fourier transform of the complex function F(z). This function and thus the structure 
of the surface layer can be reconstructed directly by the use of the inverse Fourier transform. However, it is 
necessary to know both the amplitude and phase of R{9). 

The amplitude of R{9) can be obtained from the angular dependence of the reflectivity curve measured 
in a big angular range far from the exact Bragg condition according to equation 

p R (e) = c p \R(0)\ 2 , (no) 

where C p is a polarization factor equal to C p = (1 + C 3 )/(l + C) in a double-crystal experimental scheme 
and parameter C is defined in (^0|). 

The phase a(9) of the scattering amplitude R(9) can be determined by measuring SR yield with a 
small depth of yield L y i satisfying the condition: L y i « Ld in a big angular range far from exact Bragg 
condition. This can be, for example, photoelectron yield coming from the submicron depths of the sample 
or the fluorescence yield of the foreign atoms absorbed or implanted on the top of the sample. This angular 
dependence can be obtained directly from the general equation (|52|) and is equal to 

k(A0) = l + C p \R(A9)\ 2 + 2C P \R(A9)\ |e^| /(0) cos [$ R (A6) + (p(0) + arg(e^)] , (111) 
where /(0) and ip(0) are the values of the static Debye-Waller factor f(z) and phase (fi(z) at z = 0. As 



it follows from (|ll0j) and (111) by simultaneous measurement of the reflectivity Pr{9) and SR yield k(6) 
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both the amplitude \R(A9)\ and the phase a(0) of the scattering amplitude R(9) is obtained uniquely in the 
whole angular range. The values of /(0) and <p(0) can be found from the fitting of the central part of curves 
Pr(0) and k(6), where conditions of the strong dynamical scattering are satisfied. Finally, using inverse 
Fourier transform, the structure of the damaged layer can be obtained directly. 

In the paper [ 105 1 this approach was applied for the determination of the depth profile of distortions of 
the surface layers of silicon crystals after implantation. Samples were implanted by B + ions (with energy 
of 100 keV, dose of 10 15 cm -2 ) and then annealed (T = 800°C, t — 10 min) in an oxidized atmosphere. In 
experiment the (400) reflection and CuK a radiation were used. The photoclectron emission and the x-ray 
reflectivity from the samples were measured both at large deviations from the Bragg angle (Fig. 42a) and 
in the region of total reflection (Fig. 42b). Electrons from K— shell with energies from 5.9 keV to 6.4 keV 
were excited and detected by a gas proportional counter with a resolution sw 15 -j- 20%. According to Eq. 
( |56|) an estimated depth of yield of this electrons was equal to L y i = 0.1 fim. 

It is important to note here, that reflectivity curves obtained in a standard double-crystal scheme of mea- 
surements are the superposition of coherent and incoherent scattering channels and therefore are inapplicable 
for our theoretical analysis, where only coherent part of the scattered radiation has to be used for the phase 
retrieval. This problem can be solved by measuring the "tails" of the reflectivity curves in a triple-crystal 
scheme. This arrangement makes it possible to separate the coherent part of scattering from its incoherent 
part (Tog) . 

Photocmission curve in the region of the strong dynamical diffraction (Fig. 42b) is similar to the inverted 
curve from an ideal crystal. It coincides with theoretical one with the values of the phase <p(0) = 2.3 and 
static Debye- Waller factor /(0) = 0.9. Described above algorithm for the phase reconstruction was applied 
to obtain the structural information about the surface layer of the sample from the "tails" of experimentally 
measured curves (Fig. 42a). Profiles of the deformation Ad(z)/d and Debye- Waller factor f(z) obtained 
according to this approach are presented on Fig. 43. From this Figure it is well seen, that the deformation 
in the damaged layer is positive and reaches its maximum at a depth L ss 0.3 \xm. This result correlates well 
both with the distribution of B + ions in silicon lattice under the conditions of implantation and the theory 
for defect formation in these crystals [ 107 1 . We have to note here, that all calculations were carried out in 
the limit L y i — > 0. The final depth of yield of the photoelectrons can be also taken into account, however, 
this does not change result essentially. 

Finally, in the paper [105] it was presented a method for a direct solution of the phase problem while the 
x-ray scattering in the deformed crystal. From simultaneous measurements of the photoelectron emission 
and reflectivity curves in a big angular range near exact Bragg angle the complex scattering amplitude can 
be determined completely and thereby the structure of the deformed surface layer can be reconstructed 
uniquely. 
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10 Conclusions 



Traditionally, since the discovery of x-rays, x-ray diffraction methods were used to obtain structural informa- 
tion from different type of crystals, including protein crystals [ 1 00 1 . In this review we have shown how during 
last 20 years combination of the traditional x-ray diffraction and x-ray standing wave technique has become 
one of the most sensitive methods of studying structure of real crystals. Here we have treated in detail the 
problems involved in theoretical studying and practical application of x-ray standing wave method for the 
structural analysis of the surface layers. Starting from the perfect crystal, crystals with amorphous surface 
layer were analysed. Then theory and different applications of the simplest model of deformed surface layer, 
so-called bicrystal model were presented. This model was generalized for the case of a multilayer model. 
Analytical approach for the secondary radiation yield from such crystal was analysed in details and different 
applications of this approach were discussed. It was shown, that both geometries of diffraction Bragg and 
Laue can be studied in the frame of one theoretical approach. However x-ray standing waves and secondary 
radiation yield in Laue geometry has its own peculiarities different from Bragg geometry and they were 
discussed separately. This theory was generalized for the case of the crystals with an extended deformation 
field, such as crystals with a uniform strain gradient (bent crystals) and vibrating crystals. We have also 
discussed the possibilities of the direct measurements of the phase of the scattering amplitude from deformed 
crystal and presented the way of the unique reconstruction of the structure of this deformed layer. 

During last decade new materials with high crystal quality has appeared such as High-T c superconductors, 
quasicrystals, isotopically pure materials. This opens the possibility to use XSW analysis in such crystals. 
For example, High-T c superconducting crystals Ndi^Ceo.ibCuO^s and ultrathin films of YBa2Cu307 on 
SrTiO3(001) were studied with XSW method in f54| |lOSj] . Fluorescence radiation from AlPdMn icosahedral 
quasicrystal excited by standing waves formed from x-ray reflections along a twofold symmetry axis was 
measured in a recent paper |109| , lattice constants of the isotopically pure layers of 76 Ge grown on natural 
germanium were determined while the photoelectron yield XSW studies in j78). The effect of different 
methods of growth on the interface quality of extremely thin buried Ge 5 layers in Si were studied by means 
of different diffraction techniques including XSW [ 1 10 1 . Photoelectron yield and fluorescence radiation from 
(AlAs)(GaAs) short-period superlattices was measured in 
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One interesting approach in XSW analysis was developed in ]113[ . "Kossel diffraction in perfect crystals 
can be viewed as the reverse process of x-ray standing waves with the source of radiation located inside the 
crystal. New experimental method was proposed in this paper for determination of atomic positions in the 
bulk or at a surface of a crystal by analysing the profile of angularly resolved Kossel lines. 

In the last years different theoretical approaches were elaborated for the analysis of XSW in real crystals. 
For example , x-ray standing waves in crystals with statistically distributed defects were analysed theoretically 
Analytical results for standing wave behavior in the deformed crystal with a special type of 



in [114 



Standing waves in Laue geometry for a special type of scattering with 



deformation were obtained in | 
only imaginary part of the atomic scattering factor were analysed theoretically in 116 . 

As it was shown in our review in many situations it is important to reduce the signal coming from 
the background and to enhance the surface contribution. From this point it is quite favorable to make 
experiments in grazing incidence geometry proposed for the first time by W.C. Marra, P. Eisenberger and 
A.Y. Cho 1 117 1 a nd later d eveloped as a useful tool for surface sensitive XSW measurements in a number of 
papers |l8t |ll9j |l20| , |m) (see also a book jT|]). 

In this review paper we were discussing in details formation of XSW field in the case of two-beam dy- 
namical diffraction in real crystals. This approach can be generalized for the case of multibeam dynamical 
diffraction. This was done for the first time theoretically in |}8|, 122]. Angular dependencies of the photo- 
electron yield in the multibeam region for the case of four-beam (220/400/220) and three-beam (444/335) 
dynamical diffraction in Si perfect crystals were calculated. If in the case of two-beam diffraction in one 
measurement information about atomic distribution in the sublattice can be obtained only along h-vector, 
in the case of the multibeam diffraction when second reciprocal vector g is involved in diffraction process one 
can obtain structural information simultaneously in both directions. For the first time fluorescence measure- 
ment in the case of multibeam diffraction for a perfect crystal was performed in 1 123 [ , later photoemission 
yield from a perfect Si crystal for the case of (111/220) three-beam diffraction was measured in [124, |125fl 
(see also for review of this results 126 ). Here we would like to mention also another method of direct 
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(112) 



measurement of the phase of the structure factor or total deformation field in thin surface layers without 
measuring any SR yield. This method is based on the measurement of the "weak" three-beam reflection far 
from the multibeam region simultaneously with the strong one. In this case intensity of this reflection can 
be presented in the form 

% XgO 

where a g is the angular deviation from exact multibeam Bragg region. This expression is similar to ( |77j ) 
obtained for the SR yield with the small depth of yield. This approach was for the first time applied for 
the measurement of the phase of the structure factor in perfect crystals in [ 101 1 and then generalized for the 
case of deformed crystals theoretically [122, 127 and experimentally [12S|. 

Approach described in this review paper to the calculation and analysis of the secondary radiation yield 
while the dynamical diffraction of x-rays was quite general. However we have restricted ourself mainly with 
the cases of fluorescence and photoelectron yield. Actually the potentialities of applications of x-ray standing 
waves are much broader. They include for example the study of such inelastic processes as thermal diffuse 
and Compton scattering and internal photoeffect under conditions of dynamical diffraction of the incident 
radiation. Every process has its own peculiarities and contains information of a fundamentally new type. 

For example, a unique information about high-order total photoelectron cross sections, radial wave func- 
tions and electron wave phase shifts can be obtained from XSW angular resolved and integrated photoelectron 
yield (see for theoretical and first experimental results p8l E9, 129, 130, 131]). We want to mention here 



also recent results [132] where the photoelectron yield in the field of XSW was measured from the valence 
band of the Ge and GaAs crystals giving the possibility to measure the bond polarity parameter. 

Thermal diffuse and Compton scattering in the conditions of dynamical diffraction were studied experi- 
mentally starting from the 60-th [H], [l^] (see also theoretical work Q[). In jl6| it was shown the possibility 
to determine the values of the phonon eigenvectors directly from XSW experiment. In the theoretical papers 
p5| fill unified Green function approach to the problem of the yield of the thermal diffuse scattering was elab- 
orated and it was shown that surface contribution in a special diffract ion conditions can become important . 



Compton scattering in the XSW fiel d wa s studied intensively in 80-th [133, 134, 135, 136, 137 , 138 , 139 
(see also for review of this results 



141]). In this papers it was shown, that investigation of the coherent 



Com pton scattering opens the possibility for determination of non-diagonal elements of th e density mat rix 
[136 1 and contribution of the valence electrons to the atomic factor of the x-ray scattering 1 13S , 13£, 14C]. 

Measurements of the photoelectric current in a semiconductor crystals while the dynamical scattering of 
x-rays |l7], [l^] opens the possibility to study b oth t he s tructural and the electrophysical characteristics of the 
crystal simultaneously. In a series of papers [142, 143] it was demonstrated, that the depth of yield in this 
case is determined by the diffusion length Ld of minority carriers. X-ray luminescence studies of phosphor 
crystals excited by x-ray standing waves were carried out in [144]. 

In this review paper we were mainly discussing the possibilities of generating the XSW in real crystals not 
in extreme conditions. However, there are exist some ways for generating x-ray standing waves in so-called 
extreme conditions with Bragg diffraction close to 9b = 90° and total external reflection (so-called 0— th 
order Bragg reflection). Due to extreme broadening of the Darwin curve in the first case (the half- width of 
the curve A61/2 in this case is proportional tOy/xh) this conditions of generating XSW can be effectively 
applied for less perfect for e.g. mosaic crystals (see for modifications of the dynamical theory [145, 146] and 
for XSW measurements [147, 14S, 14E, 15C] and also review paper p3|). 

The case of the total external reflection can be treated as a special case of dynamical diffraction with 
0— th order Bragg reflection. The amplitude ratio for the reflected amplitude En to the incident amplitude 
Eq can be obtained from the Fresnel formula ]151| or directly from ( p3| , |3l| ) with 9b = 0, \h = Xo an d 
approximation sin 9 w 9 



Er 
E 



Vo 2 



Xo 



Xo 



(113) 



The peculiarity of XSW at the total external reflection is the strong variation of the wavefield spacing from 
the reflection angle from about 50 A for 9 « 9 C (where 9 C = ^/\xo\ is the critical angle) to infinity for 9 — > 0. 
For experimental and theoretical results of XSW measurements in the total reflection conditions see for e. 

g. m, EM EM EM' [Oil- 
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Another way to generate an x-ray interference field with long-period wavefield spacing is the use of Bragg 
diffraction from Layered Synthetic Microstructures (LSM) or organic Langmur-Blodgett (LB) multilayers. 
Naturally, for generating standing waves these structures has to be periodic. Alternating amorphous layers 
of high— Z \ow—Z materials (e.g. W/C, Rh/C or Ni/G and etc.) provides high maximum reflectivity. The 
lattice plane spacing in such materials can typically range from 15 to 200 A. After the first demonstration 
by Barbee and Warburton [ 157 1 of the possibility to use these structures in XSW analysis it was applied 



effectively in a number of papers [158, 159, 160, 161, 162, 163 



In conclusion we can see that XSW method due to its unique properties of combination of the x-ray 
dynamical diffraction with spectroscopic technique has become a powerful tool in modern crystallography 
and even more general in condensed matter physics. 
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11 Appendix 



Secondary radiation yield from a multilayer crystal 
(analytical approach) 



In Appendix we will give the details of the analytical approach (following p3| ) for the angular dependence 
of the SR yield for a multilayer crystal considered in Section VII. Here both geometries (Bragg and Laue) 
will be discussed simultaneously. Evidently, bicrystal model considered in Section V is a particular case of 
this more general approach. 

Following this approach deformed layer with an arbitrary shape is considered as a sum of several sublayers 
(see Fig. 24) with constant parameters (deformation Ad/d and amorphization e~ w ). In this case TT 
equations in each sublayer can be solved analytically. Substituting this solutions in the general equation ( |52] ) 
angular dependence for the SR yield in analytical form for a multilayer crystal can be obtained. 

Solution of Eq. (0) with constant coefficients for the amplitude R(z, 9) in the n-th layer z n < z < z n + d n 
has the following form 



R(z,( 



F u {z, 9) _ xi - x 2 x 3 e 
F d (z,9) ~ 1 - x 3 e^ L 



(114) 



where the upper sign belong to Bragg geometry and the lower one to Laue. In this equation the following 
notations are used 



F u (z,6) =x x - x 2 x 3 e 



,F d (z,( 



1 - x 3 e 



(115) 



and 



•? - 3 



X 3 



1 

xi - 



b- 



X2 - Rd, 

xi - R„ 

X 2 - Rn 



_ e <jd n (Bragg case), 
(Laue case) 



1 



(116) 



with 



-y(0) - iya + y^^a 



2i 

Lex 



b 2 T Cl A; 



(117) 



and parameters y(9), yo, y^> and C\ arc defined in ( |25f]27| ). For the square root in Eqs. (116, 117) the branch 
with the positive imaginary part is chosen. The amplitudes R n and Rd n in (116) are the values of R{z,9) 
on the top (Az„ = 0) and on the bottom (Az n = d n ) of the n-th layer boundaries. 

Angular parameters y(9),yo,y v defined in ( p5| , ^6| ) can be written in a more general and invariant form 
through the extinction length L ex (|28|). Not taking into account the shift of the origin on the axis A9 (due 
to refraction) and defining by 69 the angular shift due to deformation in the layer we have 



y{9) = C yt A9, y = ± 



470 



-(l±/3), y v = C 



ytC- 



(118) 



In (118) parameter C y t gives connection between the angles measured in angular units (arcsec) and in 
dimcnsionless parameter y 

C yt = ng^— sin 26 B , .g = 0.4848 ■ 10 -5 . 
A 70 

If extinction length L ex is defined with the parameters of the substrate, then parameters C y t and L ex in all 
layers are the same, however, such parameters as 69, e 
from layer to layer. 



-W 



/lo and p (sec Eq. (|2?])) can, in principle, differ 
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From Eq. ( p2| ) formal solution for the direct beam Eq(z, 9) in the layer z n < z < z n + d n can be obtain 
in the following form 



E (z,6) = E (z n ,9)exp 



A70 



z 

Az n - i^- [ dz'R(z',t 

L ex J 



(119) 



Integral in the exponent with the amplitude (114) can be taken analytically giving for the direct field 
amplitude 

E Q (z,e)=T(z,9)E (z n ,e). 



Here T(z, 9) is the transmission amplitude through the layer 

T(z,9) = exp W )AW2] ( , m = ^ _ ^1. 

V 1-^3 / A70 L ex 

Intensity transfer through the layer is given by the function 



\T„ 



|T(d„,l 



where 



-M(0) 



M{9) = Im[<S>{9)] = J^(l T (3) Ta 
270 



\F d (0, 



(120) 



(121) 



(122) 



(123) 



According to (J52|) the normalized SR yield from the n-th layer in general case of multiatomic unit cell for 
the sort of atoms a and with exponential form of the yield function ( |55| ) from the entrance surface can be 
presented in the following form 



z„+d„ 



<(0) = X% / dzP yi (z)\T(z,9)\ 2 x 



1 



\E h {z,l 



2CRe 



\E (z,9)\z 

Taking into account our definitions this can be written as 



Eh\Z, 9) gJ V °(g)-W 

E (z,9f^ 



(124) 



<(0) = Xio dzP yi (z)\T(z,9)\ 2 {l + P\Y\ 2 \R(z, 



2y//3CRe \YR(z,9)e%t 



,i( V a (z)- V (z))-W a 



(125) 



In the following it will be assumed, that ip a (z) — f(z) = A(p a — const. This condition mean in fact, that 
due to macroscopical displacements the unit cell is shifted as a wh ole. S ubst ituting now previously obta ined 
solutions for the reflected R(z,9) jlTil) and transmitted T(z,9) (^2lj ^22|) amplitudes into Eq. (|l2^ ) SR 
yield from the n-th layer can be obtained 



<(0) = 



Xi 



\F d (0,i 



dze -(M{e)+„ yi )z x 



+K r \F u (z, 9)\ 2 + Re [KiF^z, 9)F u {z, 9)]} 
where functions F u (z,9) and Fd(z,9) are defined in |l!5| and 

K r = [3\Y\ 2 , Ki = 2^f(3CYF«e 1 ^. 



(126) 



(127) 
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In this equation two important parameters of the theory i.e. coherent fraction (P c a ) and coherent position 
(P") for the atoms of the sort a in the n-th layer are introduced. They are defined as following 



F c a = \e^\e~ wa , ^ = 2irP c a = A<p a + arg^), (128) 

Integral in the expression for can be taken analytically and finally for the angular dependence of the 

SR yield from the n-th layer we have 

<V) = u^^l I^i + A ^ ~ M^)] , (129) 

| -l x 3 \ 



where 



A! = l + K r \ Xl \ 2 + Re(K iXl ), 

A 2 = \x 3 \ 2 [l + K r \x 2 \ 2 + Re(K lX2 )] , (130) 
A 3 = x 3 [2(1 + K r x\x 2 ) + KiX 2 + K*x\] 



and 



i>i = (1 - e"*)//i, <k = (l- e-f*)/f 2 , ^3 = (1 - e~^)/f 3 , 

h = (M(6) + ti yi )d n , f 2 = fi± 2a r d, f 3 = fx± ad. (131) 
The total yield of the SR from the whole multilayer crystal is the sum of the yields from every layer 

N 

< ot (0) = 5>„(0)<(0), (132) 

n=l 

where z x = 1 and z n (6) = z„_i(0)|T(d„, 6)\ 2 e-^' d " . In the Bragg case the reflectivity is equal to 

P R (9) = \YR(0,6)\ 2 (133) 

and calculation is performed layer by layer from down to the top starting from the iV-th layer up to the first. 
In the Laue case transmission coefficient can be calculated according to the relationship 

N 

P T {6) = z N+1 (0) exp(^ ^dn) (134) 

and reflectivity in the same geometry is equal to 

P R (0) = \YR N+1 (6)\ 2 P T (6). (135) 

Finally, the whole algorithm for calculation of the angular dependence of the SR yield from the multilayer 
crystal in both (Bragg and Laue) geometries is defined. 
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Table 1. Values of ipo, Wo and Ad/d obtained as a result of fitting of experimental data presented on 
Fig. 16. Calculations were made for different concentrations of Ge (n~f e ) in the epitaxial film in the frame 
of bicrystal model. 



n^ e (10 19 atoms cm 3 ) 


Wo 




{Ad/d) x 10 b 





0.51(20) 


-0.84(20) 


-2.2(5) 


3.7 


0.25(20) 


-0.55(20) 


-1.4(5) 


7 


0.8(2) 


-0.09(20) 


-0.2(5) 


10 


0.8(2) 


0.37(20) 


0.9(5) 


15 


0.79(20) 


1.36(20) 


3.6(5) 
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Table 2. Experimental parameters measured by photoluminescence (PL), x-ray rocking curve anal- 
ysis (XRC) and x-ray standing wave method (XSW) for 10 nm CaF2/Si(lU) layers grown at different 
temperatures. 



Ti (°C) 


100 


550 


660 


770 


Aa/a (%) 


0.8 


1.7 


1.9 


2.2 


AE (cm" 1 ) 


-38 


-30 




100 


SX (nm) 


0.7 


1.4 




1.6 


e^(300) (%) 


-0.60 


-0.55 




0.54 


eL nu (%) 


-0.58 


-0.42 


0.007 


0.55 


A^in (arcscc) 


-640 


-550 


-320 


-25 


-W 

e 


0.80 


0.80 


0.65 


0.65 


Pxsw (A) 


2.91 


2.86 


2.83 


2.91 
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Table 3. Parameters of the layers obtained as a result of theoretical fitting of the fluorescence yield from 
Si samples implanted by Fe ions after annealing (Fig. 33b). 



No. 


a \pjm) 


Aa/rf (x 10 j 


-w 
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Fc 


TZ>b c 
*c 
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0.015 
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U.U1 
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1 A 

14 
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0.01 














24 
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0.01 


-9.6(0.09) 


0.15 


0.27(0.02) 





38 
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0.01 


-9.6(0.09) 


0.34 


0.27(0.02) 





86 
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0.01 


-9.6(0.09) 


0.41 


0.27(0.02) 





100 


7 


0.01 


-9.6(0.09) 


0.54 


0.27(0.02) 





86 
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0.01 


-9.6(0.09) 


0.64 


0.27(0.02) 





40 


9 


0.01 


-9.6(0.09) 


0.77 


0.27(0.02) 





30 


10 


0.01 


-9.6(0.09) 


0.88 


0.27(0.02) 





14 


11 


0.02 


-9.6(0.09) 


0.97 


0.27(0.02) 
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Figure 1 : Interference x-ray standing wave field as a result of coherent superposition of the incident Eo (r) 
and diffracted E/ l (r) plane waves with the periodicity d — 2n/\h\. 
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Figure 2: Crystal of thickness L, with the deformed surface layer of thickness Ld- While the Bragg diffraction 
x-rays penetrate into the crystal on the depths of the order of L ex . The depth of the SR yield is determined 
by the parameter L V i. It is shown the situation, when the condition L y i << Ld << L ex is satished. 
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Figure 3: Shcmatical view of dynamical diffraction and photoclcctron yield measurements excited by two 
coherent electromagnetic waves Eo(r), E^(r) characterised by the wave vectors ko, and polarization 
vectors eo a , and e/j CT , . The angle 29 is the angle between the wave vectors ko and k^ . The direction 
of the escaping photoelectron is determined by the vector k p .The relationship of the vectors k , k^ and h 
is indicated. 
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Figure 4: Angular dependences of the intensity of the wave held at the atomic planes (I, solid line), the 
reflectivity (Pr, dotted line), the interference absorption coefficient (fan dot-dash line), and the phase a 
(dash-double-dot line). 
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Figure 5: The values of the zero Fourier components of the susceptibilities \0r and xoi calculated for Si 
and Ge as a functions of energy E 
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Figure 6: Calculated cross-sections for photon scattering from carbon showing the contributions of pho- 
toelectric, elastic (Raylcigh), inelastic (Compton), and pair-production cross-sections. Also shown are the 
experimental data (open circles). From Rcf Even for the light clement carbon, the photoelectric 

cross-section is clearly dominating in the x-ray energy range up to more than 100 keV. 
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Figure 7: Part of x-ray radiation while scattering on crystal is absorbed in a layer with the thickness dz at 
the depth z and gives rise to the yield of the secondary radiation (photoclcctron e~ , fluorescence hv, etc.). 
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Figure 8: Reflectivity Pr(u) (a) and photo-emission K(y) (b) curves calculated for a Si crystal with a de- 
formed surface layer under conditions of (111) diffraction of a plane monochromatic wave of CuK a radiation. 
Also shown (c) the profile of deformation Ad/d(z) in the surface layer. From Ref. jfidjj. 
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Figure 9: Left: Angular dependence of the yield of GeK a fluorescence and of x-ray reflection in (111) 
diffraction of MoK a radiation. The fluorescence data are measured with a Si(Li) detector situated at 
a glancing angle to the surface, the value of which is indicated in milliradians. From Ref. J5^/; Right: 
The angular variation of the perfect Si (111) reflectivity and photoclectron yields while diffraction of x-ray 
radiation with energy E = 15.1 keV. Photoclectron curves correspond to electrons with different energy 
losses AE (AE = (A), 2.5 (B) and 5 keV (C) respectively). From Ref. j^Bj. 
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Figure 10: XSW measurement of CuK fluorescence radiation from a Cu crystal. Measured (+++) and 
calculated (■ ■ ■) Cu(lll) reflection curve; (•) experimentally determined fluorescence yield; ( — ) fit to the 
experimental points. Two different dislocation densities DD and takeoff angles a were used. Fluorescence 
yield curve for the big exit angle (a — 10°) coincide with the reverse reflectivity curve. From Rcf. M. 
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Figure 11: Angular dependences of the secondary radiation yield calculated for Si (400) reflection and 
CuK a radiation for different values of the ratio L ex /L yi equal to: 12 (1), 7.2 (2), 3.6 (3), 1.8 (4), 0.9 (5), 
0.45 (6), and 0.225 (7), in a system with a symmetric monochromator. From Ref. JIT 
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Figure 12: Top: Schematic view showing the position of the noncentrosymmetric GaAs(lll) diffraction 
planes (dashed lines) relative to the Ga atoms (open circles) and As atoms (closed circles). Parameter P c 
locates positions in this structure in the (111) direction relative to a (111) diffraction plane, which is shifted 
by an amount A relative to a ccntrosymmetric (111) diffraction plane; Bottom: Angular dependence of the 
wave field intensity in a GaAs(\\l) crystal at the Ga (Ig<i) &nd As (Ias) atomic sites for x-ray diffraction 
with an energy of 10.372 keV (solid line) and 11.868 keV (dotted line) for a— polarization. The energies of 
radiation are chosen to be 5 eV above the GaK and AsK edges. From Ref. ffij . 
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Figure 13: AnguJar dependences of the photoelectron yield for Ge(220) diffraction, CuK a radiation with 
different thickness of amorphous layers on the surface: (1), 50 (2), 100 (3) 250 (4) 450 (5) nm. From Rcf. 



71 




0% ■■ - i— b i j ■ ■ i i 

REFLECTION ANGLE 

Figure 14: The angular variation of the 0.6-fim S1O2/ Si(\l\) photoelectron yields for energy regions A-D 
of the emission spectrum (shown in the inset). The solid lines are theoretical curves. From Ref. |53/. 
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Figure 15: The profile of deformation (Ad(z)/d) in a crystal for the bicrystal model. Intcrplaner distance 
in the layer d + Ad differs from that in the layer on a constant value Ad. In this model interface is a sharp 
(step) function located on the depth Ld- 
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Figure 16: Photocmission angular dependence for Si single crystals with autoepitaxial Si Rim, doped with 
B and Ge with varying Ge concentrations: (curve a), 3.7 x 10 19 (b), 7 x 10 19 (c), 10 20 (d), 1.5 x 10 20 (e). 
From Ret. 
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Figure 17: Angular dependence of InL a fluorescence yield from the film of Ino.5Gao.5P (thickness of the 
film ~ 100 A) grown on the surface of GaAs (111) single crystal (upper curve) and reflectivity (bottom 
curve). Points are experimental data and lines result of the theoretical fitting. From Rcf. 
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Figure 18: Reflectivity (Top) and angular dependence of InL a (upper curve) and PK a (bottom curve) 
fluorescence yield (Bottom) from the film of Ino.5Gao.5P (thickness of the film ~ 0.6 urn) grown on the 
surface of GaAs (111) single crystal. The shadowed region on the reflectivity curve is the region, where 
angular dependence of the fluorescence yield was measured. Points are experimental data and lines result of 
the theoretical fitting. From Ref. 
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Figure 19: (Left:) The (222) x-ray reflectivity curves of Cai 7 ^ epitaxial 10 nm layers grown at different 
temperatures; (Right:) Angular dependence of the total photoclcctron yield k (upper curve) and reflectivity 
Pr (bottom curve) for the (888) diffraction of CuK a radiation in a perfect GGG garnet single crystal. Points 
are experimental data and lines are theoretical fitting. From Ref. fffdj . 
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Figure 20: (a) Schematic arrangement of the planes consisting from the different sort of atoms in positions 
a, c and d in a single garnet crystal in the direction [111] in the unit cell; (b) Angular dependence of the total 
photoclectron yield k (upper curve) and reflectivity Pr (bottom curve) for the (888) diffraction of CuK a 
radiation in a perfect GGG garnet single crystal. Points are experimental data and lines are theoretical 
fitting. From Rcf. fiFJj. 
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Figure 21: Experimental results of the photoelectron yield k (upper curve) and reflectivity Pr (bottom 
curve) for (888) diffraction of CuK a radiation in a perfect GGG garnet crystal with an epitaxial film of FYG 
garnet crystal with the thickness 2 urn. From Rcf. f^/ . 
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Figure 22: (Left:) Angular dependence of the fluorescence yield of NdL a (upper curve) and GaK a (middle 
curve) for the (444) diffraction of MoK a radiation (reflectivity is shown on the bottom curve) for the perfect 
GNG garnet crystal; (Right:) Angular dependence of the YK a fluorescence yield from the Rim (upper curve) 
and GaK a fluorescence yield from the substrate (middle curve) for the (444) diffraction of AgK a radiation 
(reflectivity is shown on the bottom curve) for the perfect GGG garnet crystal with an epitaxial him of FYG 
crystal with the thickness 1.6 fim. Points are experimental data and lines are theoretical htting. From Rcf. 
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Figure 23: Projection of the dodccahcdral positions c z (full circles) and c xy (open circles) on the plane 
(001). The (004) x-ray standing wave is also shown schematically. From Ref. pdjj. 
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Figure 24: Calculated angular dependencies of the reflectivity Pr{&) (a) and fluorescence yield of the Bi 
(b) for the (004) diffraction of MoK a radiation on Y^- x Bi x F 'e^O^ thin film for the values of p from to 1 
in the angular range of the x-ray diffraction on the him. From Rcf. |7q/ . 
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Figure 25: Experimental (points) and calculated (line) diffraction curves in the angular range of the x-ray 
diffraction on the film and on the substrate. From Ref. f7g|/ . 
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Figure 26: Experimental (points) and calculated (lines) for the value of p = 0.44 angular dependencies of 
the fluorescence yield of Bi 3+ (a), Y 3+ (b), Fe 3+ (c) ions from the film and Gd 3+ ions from the substrate (d) 
and the x-ray reflectivity (e). For comparison the fluorescence yield curve for p — 0.33 (uniform distribution 
of Bi 3+ ions over c z and c xy sites) is also shown in (a). From Rcf |7fi|/ . 
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Figure 27: (Top:) Experimental set up for SR yield measurement in Laue geometry: C\ crystal monochro- 
mator, C2 crystal under investigation, Nal detector measuring reflectivity, Si(Li) detector measuring fluo- 
rescence yield. (Bottom:) Distribution of x-ray standing wave field in a crystal for Laue geometry. Nodes 
of the weakly absorbing Geld (upper curve) coincide with the atomic planes and maximum of Geld intensity 
for highly absorbing (bottom curve) coincide with atomic planes giving maximum absorption of this Geld on 
the thickness of the crystal L. 
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Figure 28: Angular dependence of the GeK a fluorescence yield from the Si crystal uniformly doped with 
Ge with different thickness of the crystal. Points are experimental data, lines are calculations. (Left:) Case 
of the thick crystal (L = 2.2 mm). Solid line is the calculation for substitutional impurity atoms (coherent 
position P c m = 0); dashed curve is the calculation for impurities with coherent position P,, 111 = 0.15. The 
x-ray reflectivity curve Pr is also shown. (Right:) Case of the thin crystal (L = 0.49 mm). Curve 1 is the 
calculation for the substitutional impurity atoms and curve 2 is the calculations for the randomly distributed 
impurities. From Rcf. $$3jj - 
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Figure 29: Theoretical dependencies of the maximum fluorescence yield for substitutional impurity atoms 
in silicon crystals with different thickness. Solid lines are 111 reflection and dashed lines are 220 reflections. 
FromRcf. J§|/. 
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Figure 30: Calculated angular dependencies of the impurity fluorescence yield from a thin Si crystal for a 
different impurity positions. From Rcf. 
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Figure 31: Angular dependence of the GeK a fluorescence yield from the Si crystals with an epitaxial layer 
with thickness 0.35 mm. (Left) Angular region corresponding to the x-ray diffraction on the substrate. 
Curve 1 is the calculation for the substitutional impurity atoms; curve 2 is the calculation for the randomly 
distributed impurities; circles are the experimental data. (Right) Angular region corresponding to the x- 
ray diffraction on the layer. Fluorescence yield is presented on the upper curve and reflectivity on the 
bottom. Dashed lines are calculations for substitutional impurity atoms (1 ) and for the randomly distributed 
impurities (2). From Rcf. Jgj|/ . 
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Figure 32: Schematic view of the SR yield from the multilayer crystal while the Bragg diffraction of x-rays. 
The n-th layer of such crystal is characterized by constant values of deformation Ad/d and amorphization 
e~ w . The x-ray scattering from each layer is given by its reflection R n and transmission T n amplitudes. The 
SR yield from the n-th layer and the atoms of the sort a is determined by the values of coherent fraction F c a 
and coherent position P". 
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Figure 33: The fluorescence FeK a yield for silicon crystals implanted with Fe ions: just after implantation 
(a) and after annealing at 750° (b). Points are experimental data and lines result of theoretical fitting (fitting 
parameters are summarized in Table 3). The reflectivity curve (c) for the sample (Si (111 ) reflection, CuK a 
radiation ) is the same before and after annealing. From Rcf. pdJ. 
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Figure 34: (a) Damage density in silicon due to 80 keV Fe ion implantation for the sample just after 
implantation(circles) and after annealing (triangles), (b) Profile of Fe impurity concentration for the sample 
just after implantation, (c) Profile of Fe impurity concentration for the sample after annealing. From Ref. 
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Figure 35: Schematic view of the multilayer model used for the calculation of the x-ray waveheld in a crystal 
with USG. Thickness of each layer is Az = L ex /(mC). The angular deviation parameter yi(9) is constant in 
each layer and differs from layer to layer by the magnitude Ay(8) = l/m. The amplitudes Eq, R n , T n are 
the amplitudes of incoming, reflected and transmitted beams from the n-th layer respectively. 
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Figure 36: Angular dependences of the reflectivity, fluorescence and photocmission yield calculated for a 
cylindrically bent Si (400) crystal, CuK a radiation, at different values of the curvature parameter C: (a) 
(perfect crystal), (b) 0.1, (c) 1.0, (d) 10. From Ref. 
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Figure 38: The two-dimensional distribution of the standing wave intensity I(z,y) calculated at depths 
z = -j- 2L ex for the curvature parameter C = 1 . In the inset the angular dependence of the intensity for 
z = is plotted. (Top) Intensity values calculated on diffraction planes (Ay>° = 0). (Bottom) Intensity 
values calculated between diffraction planes (Aip^ — n). From Rcf ^dj . 
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Figure 39: The angular dependence of the x-ray reflectivity (a) and the x-ray standing waves (b and c) 
calculated on the surface of the Si (400) crystal without deformation (left figure) and with a USG (right 
figure). The XSW curves correspond to different position (the phase shift Aip®) of the impurity atoms on 
the surface of a crystal. From Rcf. j9i|/. 
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Figure 40: Schematic view of the secondary radiation yield from a vibrating crystal. This crystal is presented 
as a set of compressed and decompressed layers. 
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Figure 41: The angular dependence of the photoelectron yield and reflectivity calculated for the values of 
the amplitude hw equal to 0.1 (curve 1), 0.5 (curve 2), 0.9 (curve 3). Figures (a) and (c) correspond to 
satellites with N = =pl and (h) to the main reflex (N = 0).From Ref. |9j| /. 
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Figure 42: Photoemission yield (upper curves) and x-ray reflectivity (bottom curves) measured at big 
deviations from exact Bragg condition (Top) and in the region of the strong Bragg diffraction (Bottom). 
Points are experimental data and lines are theoretical calculations. Curve 3 on Bottom Figure correspond 
to the photoemission yield from an ideal silicon crystal. From Ref. Ji05f. 
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Figure 43: Profiles of the deformation Ad/d(z) and of the static Debye-Waller factor f(z) = exp(— W) 
obtained as a result of the phase retrieval from the angular dependence of the photoemission yield and 
reflectivity (Fig. 42a). From Re f. JlQgj/. 
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